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Abstract 

Quantum electrodynamics is the first successful and still the most successful quan- 
tum field theory. Simple atoms, being essentially QED systems, allow highly accu- 
rate theoretical predictions. Because of their simple spectra, such atoms have been 
also efficiently studied experimentally frequently offering the most precisely mea- 
sured quantities. Our review is devoted to comparison of theory and experiment in 
■ the field of precision physics of light simple atoms. In particular, we consider the 



Lamb shift in the hydrogen atom, the hyperfine structure in hydrogen, deuterium, 
helium-3 ion, muonium and positronium, as well as a number of other transitions in 
positronium. Additionally to a spectrum of unperturbed atoms, we consider annihi- 
lation decay of positronium and the g factor of bound particles in various two-body 
atoms. Special attention is paid to the uncertainty of the QED calculations due to 
the uncalculated higher-order corrections and effects of the nuclear structure. We 
also discuss applications of simple atoms to determination of several fundamental 
constants. 
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1 Introduction 



Several atoms play basic roles in modern physics and, in fact, very different 
roles. A unit of time, the second, is defined via the hyperfine interval in the 
cesium atom, while the atomic mass unit and the Avogadro number are defined 
via the mass of a carbon atom. These two atoms are significant for our system 
of units, SI. In addition, there are some favorite atomic systems where the 
basic laws of Nature find their expression in the most transparent way. These 
simple atoms, most of which consist of two bound particles, have been crucial 
for our understanding of key moments of modern physics and their study is 
still of essential interest and importance. 

The simplicity and harmony of the theory of bound systems have been tempt- 
ing and challenging for a while. Johannes Kepler believed the Solar planetary 
system to be governed by a harmony of discrete numbers via geometry, trying 
with the so-called Platonic or Regular solids. He failed to verify that kind of 
the harmony and discovered instead some regularities in the planetary orbital 
motion known now as Kepler's laws. His discovery became a milestone in the 
development of theory of gravitation and mechanics and eventually led to the 
establishment of new mechanics (classical mechanics). 

Three centuries after the Kepler's time, a planetary model was suggested for 
atoms. Meantime certain regularities in the spectrum of the hydrogen atom 
were discovered. Those regularities, like the Kepler's laws, led once again to 
the establishment of new mechanics (quantum mechanics) and simultaneously 
realized the Kepler's dream of the harmony of numbers governing the orbital 
motion. 

By now we know that a quantity describing a classical object can be of an 
arbitrary value while in the quantum case only discrete values are possible 
for some quantities. And this is how certain integer numbers enter the basic 
equations of modern physics. 

Working on a new mechanics, a new dynamics or a new model, one used to try 
first to apply it to some 'simple' objects. The simplest object is a free particle. 
However, only a limited number of properties of a free particle can be derived 
ah initio and studied with a high accuracy. Study of simple atoms opens a 
broad field for possible effects and applications. A two-body atomic system 
such as the hydrogen atom is a natural object to verify a new model or to test 
a new approach. Studies of the properties of the hydrogen atom have already 
served to establish the so-called 'old quantum mechanics' (the Bohr theory), 
modern nonrelativistic quantum mechanics, relativistic quantum mechanics 
(based on the Dirac equation) and quantum electrodynamics (QED), which 
was the first successful quantum field theory. Perhaps, we should even say that 
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QED is the only quantum field theory which is successful for a really broad 
range of problems from atomic spectra to scattering, from low energy, related 
to microwave radiation, to high energy phenomena with hard annihilation and 
bremsstrahlung, from nano- to giga- electronvolt. 

Figure 1 shows several crucial contributions to hydrogen energy levels. We note 
here that one of reasons for choosing a non-relativistic equation by Schrodinger 
over a relativistic Klein-Gordon-Fock equation was an incorrect description by 
the latter of the fine structure effects in the latter. Another remark on impor- 
tance of the hydrogen atom for QED is that the anomalous magnetic moment 
of an electron was first discovered by Rabi and his colleagues [1] as an anomaly 
in the hyperfine structure of hydrogen. Immediately that was interpreted as 
a possible anomaly related to a free electron and only afterwards was that 
confirmed by a direct experiment. A historic overview of the 'contribution' of 
the hydrogen atom to modern physics can be found in [2]. 
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Fig. 1. Some low- lying levels in the hydrogen atom (not to scale). The hydrogen 
nlj(F) levels are labelled by the values of the principal quantum number n, orbital 
moment I, electron angular momentum j and atomic angular momentum F = j + 1, 
where / is the nuclear spin. The gross structure (An ^ 0) is well explained by 
the Bohr theory (so-called 'old quantum theory') and Schrodinger theory which 
also predicts the hyperfine structure (AF ^ 0). The fine structure (Aj ^ 0) is 
explained by the Dirac theory while the Lamb shift (A/ ^ 0) is a result of bound 
state QED effects. Meanwhile effects of quantum electrodynamics for free particles 
are responsible for the a/2ir anomaly in the hyperfine structure. 

One can expect that the simplest atoms are the easiest for a theoretical pre- 
diction. That is true only in part. A simple atom usually possesses a simple 
spectrum which has relatively simple and transparent properties. The atomic 
energy levels are often perturbed during measurements by various external 
factors such as a residual magnetic or electric field. Because of the simplicity 
of the spectra, there is a good chance to understand the influence of those 
external factors on simple atoms. As examples of the highest-precision exper- 
iments on simple atoms, we remind here that the hyperfine interval in the 
ground state of hydrogen was the most accurately measured physical quantity 
for a few decades and today the Is — 2s interval in the hydrogen atom is among 
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the most precisely measured values. 



Decade after decade, theorists and experimentalists investigated simple atoms. 
As a result their theory is the most advanced atomic theory and it has to 
compete with very sophisticated experiments. The theory of simple atoms goes 
now far beyond non-relativistic quantum mechanics with the reduced mass of 
an electron. One has also to take into account relativistic effects, recoil effects, 
quantum electrodynamics, effects of the nuclear spin and nuclear structure. 

In an early time of modern physics a list of simple atoms consisted of hydrogen 
only and later deuterium, neutral helium and helium ion were added. Now the 
list is much longer and quite impressive. It also contains tritium, hydrogen- 
like ions of almost all elements (and other few-electron ions). Artificial two- 
body atoms are the easiest to produce and the list also includes muonium 
and positronium, muonic atoms, pionic, kaonic and antiprotonic atoms, exotic 
bound systems of two unstable particles (such as pionium and 7r/i-atoms) and 
antihydrogen. 

Often accuracies of theory and experiment are not compatible. However, there 
is a broad range of effects, for which theory and experiment approach the same 
high level of accuracy. The study of such effects forms a field called precision 
tests of bound state QED, which is reviewed in part here. 

A number of sources have contributed to uncertainty of such tests, and the 
current accuracy of QED calculations for free particles and two-body atoms 
is not a limiting factor for QED tests. The accuracy of the tests is limited by 
one of the three other sources: 

• an experimental uncertainty; 

• an inaccuracy of taking into account effects of the strong interactions; 

• an uncertainty due to a determination of the fundamental constants. 

The latter enters consideration because theory is not in a position to give itself 
any quantitative predictions. It provides us with some expressions containing 
values of certain fundamental constants, such as the Rydberg constant Roo, 
the fine structure constant a, the electron-to-proton mass ratio etc. To make a 
prediction, one needs first to determine the fundamental constants by extract- 
ing their values from some other experiments. Thus, theory serves as a bridge 
between different experiments. That makes the determination of fundamental 
physical constants to be another important part of precision physics of simple 
atoms. 

The contemporary situation with the QED uncertainty being below the total 
uncertainty of any QED test is a result of significant theoretical progress for 
the two last decades. Twenty or thirty years ago the QED uncertainty was 
often the dominant one. 
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The precision tests of QED form a multidisciplinary field involving atomic, 
nuclear and particle physics, laser spectroscopy, frequency metrology, acceler- 
ator physics, advanced quantum mechanics, quantum field theory etc. Those 
tests are rather not to test QED itself, but to check the overall consistency of 
the results and methods from different fields and in part to search for possible 
new physics beyond the Standard Model. In doing more and more advanced 
QED calculations, we also need to verify our approaches to deal with infrared 
and ultraviolet divergences, renormalization and bound state problem for few- 
body atoms in the case of high order terms of perturbation theory. 

As already mentioned, the simplest object to test any theory is a free par- 
ticle. A study with free leptons provides us with a possibility of testing the 
QED Lagrangian. The most accurate data for a free particle are related to the 
anomalous magnetic moments of an electron and a muon. In the former case 
the limitation of the accuracy of a comparison of theory versus experiment 
is due to a determination of the fine structure constant a while the latter 
QED test is limited by effects of strong interactions entering via hadronic in- 
termediate states for the vacuum polarization effects and by the experimental 
accuracy. The QED theory of the anomalous magnetic moment is quite ad- 
vanced including many hundreds of multiloop diagrams (up to the four-loop 
level for the electron and five- loop level for the muon). That differs from a so 
called bound state QED, a quantum electrodynamics theory of simple atoms, 
which deals mainly with relatively few one-loop and two-loop diagrams, but 
those are not for free particles, but for the particles bound by the Coulomb 
field. Three-loop contributions are rather below the uncertainty of most of the 
bound QED calculations and if necessary can be calculated neglecting binding 
effects. 

These two QED theories, the free QED and the bound state QED, are very 
different in their approaches, problems and applications and it is worth to 
consider their tests separately. The bound state problem makes all calculations 
more complicated. Intending to perform a calculation within the free QED, one 
can in principle find needed basic expressions in a number of textbooks. On 
the contrary, the bound state QED is not a well-established theory and there 
are no published common prescriptions for the relativistic quantum bound 
problem. It involves different effective approaches to solve the two-body bound 
problem. 

Precision tests of the bound state QED offer a number of different options: 

• The approaches for the bound problem can be essentially checked with the 
low-Z two-body atomic systems like hydrogen and deuterium, neutral he- 
lium and helium ions, muonium, positronium, etc. At low value of the nu- 
clear charge Z the binding energy is of order of (Za) 2 mc 2 and it is essentially 
smaller than the energy related to the rest mass mc 2 of the orbiting parti- 
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cle. That is a regime of a weak coupling in contrast to the high-Z physics. 
The weak coupling allows efficient use of a perturbation expansion over the 
Coulomb strength Za. Many contributions of higher-order in Za are calcu- 
lated diagrammatically and thus are closely related to other QED diagrams 
(e.g., for scattering theory). 

• Studies of high-Z ions are related to a strong coupling regime, however, 
it is unlikely to provide us with more information on bound state QED 
because of substantial contributions due to the nuclear structure. Such an 
investigation is rather useful for testing different nuclear models. However, 
in some particular cases, atomic systems with not too high Z can give some 
important information on higher order terms of the QED Za expansion. Im- 
portance of particular atomic systems also depends on the energy interval 
under study. For instance, the hyperfine structure (HFS) interval depends 
more strongly on the nuclear-structure effects than the Lamb shift, and 
the related calculations involve more details of the nuclear structure. The 
accuracy of the nuclear-finite-size contribution is often claimed to be very 
high, however, their uncertainty is customarily not well verified. It is often 
estimated from a variation of the result for an energy shift after applica- 
tion of certain models of the charge distribution while the charge radius is 
fixed. However, to claim such a calculation as an ab initio result, one has 
to reconsider first the original scattering and spectroscopy data on the nu- 
clear form factor and related nuclear theory and to check how strongly they 
constrain the shape of the distribution in general and its parameters and 
afterwards to question how much the nuclear-size correction varies within 
allowed shapes and parameters. That may be done only on a nucleus-by- 
nucleus basis. Lack of such a consideration in study of high-Z study reduces 
importance of the bound-state QED calculations, which are in such 

one more confirmation of applicability of the phenomenological model of 
the charge distribution. 

• Studies of few-electron atoms involve electron-electron interactions. In the 
case of high Z ions, the electron-electron interaction is proportional to a, 
while the interaction between an electron and the nucleus is proportional 
to Za. If the ion is highly charged with only few electrons remaining, the 
electron-electron interaction can be treated as a perturbation. As a result, 
at high Z the electron must be treated relativistically, i.e., no expansion over 
Za can be used, however, the treatment of the electron-electron interaction 
as a perturbation leads to the 1/Z expansion. In the case of light atoms, 
electrons are essentially nonrelativistic but the electron-electron interaction 
is compatible with the nucleus-electron interaction. The few-electron atoms 
(like, e.g., the neutral helium atom or heavy lithium-like ions) is the only 
case when the uncertainty of the QED calculations used to be a limiting 
factor for a comparison of theory versus experiment. 

• There are some other two-body atoms under investigation. They contain a 
muon or a hadron as an orbiting particle. The orbit in a muonic atom is much 
closer to the nucleus than in the case of a conventional atom with an orbiting 
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electron, and the muonic energy levels are much more sensitive to the nuclear 
structure. In the case of a hadronic atom, the orbit lies even lower than in a 
muonic atom, while the interaction of the orbiting particle and the nucleus 
is dominated by the strong interaction, in contrast to the electromagnetic 
interaction in conventional (i.e., electronic) and muonic atoms. Different 
exotic or muonic atoms offer a unique opportunity to study certain particle 
properties by spectroscopic means with high precision. 

Frequently high-Z spectroscopy is quoted as a QED test at a strong field. 
However, that is not exactly the case. A value of macroscopic meaning, such 
as the electric field strength E, has not much sense inside an atom. Other 
details are more important. For example, the strength of the field can be 
characterized by the average potential energy 

{Zafmc 2 

\ U ) ~2 > W 



which increases with the nuclear charge Z and the mass of the orbiting particle 
m and decreases with the principal quantum number n. The strongest field 
is related to high-Z atoms with an orbiting particle, heavier than electron, 
at the ground state. Muonic atoms have been studied for a while and with 
m M ~ 207m e they offer a test at a field stronger than in electronic atoms and 
at shorter distances 

(r) ~ — ^— n 2 . (2) 



The distance (r) (or a related characteristic value of the momentum transfer 
(q) ~ ft/( r )) i s another important characteristic of the electric field inside an 
atom. 

We also note that a characteristic value of the potential (U) ~ Za/(r), the 
distance (r) and the strength of the field |E| ~ (U) / (r) strongly depends 
on a particular quantity under study. What we present above in Eqs. (1) 
and (2) is related to the leading contributions to atomic energy levels. Higher 
order corrections, such as QED corrections, may involve various short- distance 
effects with characteristic distance (r) ~ h/mc and momentum transfer of 
about mc. They correspond to higher field than long-distance effects. 

The case of strong fields at short distances in which a characteristic momentum 
transfer is higher than 2m e c leads to an enhancement of polarization effects 
of the electron-positron vacuum. That makes the physics of muonic atoms 
quite different from that of conventional atoms. High-Z ions offer another op- 
portunity — a strong coupling regime with the binding energy comparable to 
mc 2 which implies a completely relativistic consideration. The strong-coupling 
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regime is very different from perturbative weak-coupling one. The calculations 
are more complicated. They may be of interest for study of methods to be ap- 
plied in the case of strong interactions. Some of the high-Z results are impor- 
tant also for high-accuracy \ow-Z physics. However, one has to remember that 
the strong coupling is rather not a problem of QED, but a problem of its ap- 
plication. In this sense, the stronger the coupling, the less 'simple' the system. 
To study pure QED problems in more detail we should prefer a weak-coupling 
situation. 

In muonic atoms the coupling constant and other parameters can take quite 
different values. While for Z — 1,2 the states of interest are n = 1,2, the 
principal quantum number n for medium Z may be higher than 1 or 2, which 
used to be common for the high-Z experiments with conventional atoms. 

However, in both cases (muonic/exotic atoms and high-Z ions), understand- 
ing the nuclear properties is needed for any theoretical predictions. For in- 
stance, let us look at studies of the Is Lamb shift in hydrogen-like uranium. 
Recently, the experimental result was improved [3]. The experimental uncer- 
tainty allows to check one-loop radiative corrections, but not two-loop effects 
calculated few years ago [4]. Those two-loop corrections are an object of in- 
tensive theoretical study and are of great interest (see Sects. 5, 7 and 13). 
The finite-nuclear-size uncertainty for the Lamb shift in U 91+ is estimated at 
a level approximately tenfold below the experimental uncertainty. However, 
the uncertainty of the result was obtained [5] (see also [6,7]) by comparison 
of two distributions of the nuclear charge, which were the homogenous spher- 
ical distribution and the Fermi distribution. The value of the mean square 
radius was fixed as (r 2 ) 1 / 2 = 5.8604(23) fm for both distributions. However, 
this result was obtained in Ref. [8] from muonic uranium spectroscopy sug- 
gesting a modified Fermi distribution, which is different from both distribu- 
tions applied in [6]. It was stated [8] that the uncertainty presented there was 
of pure statistical nature, while the model-dependence had not been studied 
and related systematic error was not included. Apparently, the characteristic 
atomic momentum in muonic uranium is much higher than in the conven- 
tional hydrogen-like uranium ion and muonic spectra are substantially more 
sensitive to the nuclear-shape effects than conventional ones. If one expects 
that a comparison of the homogenous spherical distribution and the Fermi 
distribution leads to a plausible estimation of the finite-nuclear-size effects 
(which should be verified), that should be applied first to muonic atoms to 
learn the systematic error of the mean square radius. Then, with a value of 
the radius related to each distribution, one should calculate the energy levels. 
That should substantially increase uncertainty. This example shows that how 
fragile QED tests with high-Z ions can be and how much work should be 
additionally done to develop them. 

A purpose of this paper is to give a brief review on precision physics of simple 
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atoms related to the accurate tests of quantum electrodynamics for bound 
states and the determination of fundamental constants. Because of that, we 
focus our considerations mainly on light hydrogen-like atoms (hydrogen, deu- 
terium, helium ion, muonium, positronium) and some medium Z ions, where 
the nuclear structure and hadronic effects are not too significant and allow a 
crucial test of advanced QED theory with high order contributions. 

We distinguish here QED and bound state QED, which example is an ap- 
plication of QED to the simplest atoms. Studying less simple atoms we deal 
not with just bound state QED, but with its realization for not-too-simple 
atoms. The additional problem may be related to strong field, strong cou- 
pling, crucial effects due to the nuclear structure, electron-electron interaction 
in few-electron atoms etc. Definitely, a number of investigations of less simple 
atoms are of interest. However, dealing with them is to go beyond the simplest 
bound state QED. 

We note, that the light hydrogen-like atoms are the most attractive from a 
theoretical point of view. They involve neither electron-electron interactions, 
nor strong-coupling effects or so. We consider in the next section, what is the 
most favorite choice for experimental accuracy options in testing bound state 
QED for hydrogen-like atoms. That is related to the light atoms. They are 
also favorite in principle for theoretical accuracy, being the simplest atomic 
systems. That does not mean that study of other atoms are out of interest. 
First of all, what is important is not just an atomic system, but a certain 
transition there. As one can see in this review, certain transitions, or combi- 
nations of certain quantities, related to different transitions, may offer various 
theoretical or experimental advantages. 

Because of the simplicity of simple atoms and multidisciplinary nature of the 
precision tests of the bound state QED we have tried to review the results as 
simply and briefly as possible in order to make the paper readable by non- 
experts. Detailed reference information on crucial theoretical contributions is 
collected in tables. 

While considering QED, there is always the problem of selecting units. From 
a theoretical point of view one should like to simplify equations by applying 
relativistic units in which h — c — 1, or using the atomic units. Meanwhile, 
from an experimental point of view the equations should be expressed in units 
convenient for measurements. In our paper we choose a compromise. All results 
are presented in the units of the SI. However, we present, when possible, most 
of the results using natural constants as a kind of units. For example, to define 
the so-called Fermi energy which is a result of the non-relativistic interactions 
of the electron and nuclear magnetic moments, we write for hydrogen 



and thus the proton and electron magnetic moments are explicitly expressed in 
units of the Bohr magneton. The other factors do not directly contain electron 
charge e, but only the fine structure constant a, which does not depend on 
a choice of 'macroscopic' units (in which it may be defined as a = e 2 /hc, 
a = e 2 /4irhc, a = e 2 /i.Tieohc depending on the definition of the unit for 
the electric charge). Still, here we make an exception for numerical values of 
atomic energy, which are always expressed in frequency units, i.e., in terms of 
E/h. This is because it is widely preferred to write equations for energy, while 
the actually measured quantities are the transition frequencies, spectroscopic 
linewidths and decay rates. The most frequently used notations are explained 
in Appendix A. 

More details on physics of hydrogen-like atoms can be found in: 

• various basic questions in books [9,10]; 

• an overall review with an extended comparison of theory and experiments 
related to sixties and early seventies in [11]; 

• minireviews on particular questions in [10,12,13,14]; 

• review on theory of light hydrogen-like atoms in [15,16]; 

• original results presented at Hydrogen Atom conferences and on Interna- 
tional conferences of Precision Physics of Simple Atomic Systems (PSAS) 
in [17,13,14,18,19]. 

The books [13,14] published in series Lecture Notes in Physics, volumes 570 
and 627, are also available on-line. The recent PSAS conference on simple 
atoms took place in early August 2004 in Brazil as a satellite meeting to 
the International Conference on Atomic Physics (ICAP). The coming PSAS 
meeting is scheduled for June 2006 in Venice. 

A few problems related to our paper are not presented here in detail. 

• We consider the fundamental constants only in connection to simple atoms 
and QED. More detail on fundamental constants can be found in [20,21,22,23,24]. 

• Heavy few-electron ions, which are of a great interest for study of applica- 
tion of bound state QED to strong-coupling and few-electron systems, are 
reviewed in, e.g., [6,25,26]. 

• Recent progress with exotic and muonic atoms is presented in detail in 
[27,28,29,30]. Study of such atoms are not of a big interest because of QED, 
on contrary, they deliver us a crucial information on other part of physics, 
namely, particle and nuclear physics. 

Most of this review was ready before new results on the fundamental constants 
[23] became available and through out the paper we compare the QED-related 
results on the fundamental constants with a previous set of the recommended 
constants [21]. We also note that a substantial part of QED results under 
review appeared between the publications of two recommended sets (their 
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deadlines for collecting the input data were 1998 and 2002) and most of recent 
results have been accommodated in [23]. 

We start our review with an introductory discussion of spectra of simple atoms 
and related basic QED phenomena and next consider QED tests with hydrogen 
(the Lamb shift and hyperfine structure) and other light atoms. We discuss 
studies of pure leptonic atoms such as muonium and positronium. In addition 
to spectra we consider the magnetic moments of bound particles in various 
two-body atomic systems. The concluding sections of the paper are devoted 
to fundamental constants and problems of the bound state QED. 



2 Spectrum of simple atoms and main QED phenomena 



Let us discuss the spectrum of simple two-body atoms in more detail. The gross 
structure of atomic levels in a hydrogen-like atom comes from the Schrodinger 
equation with the Coulomb potential and the result is well known: 

{Za) 2 m R c 2 

En = — 2 , (3) 



where Z is the nuclear charge in units of the proton charge, ttir is the reduced 
mass of the atomic orbiting particle (mostly, an electron) 

Mm 

™R = tt • (4) 

M + m y ' 



Here, m and M are masses of the orbiting particle and the nucleus. 
There are a number of various corrections: 

• relativistic corrections (one can find them from the Dirac equation); 

• hyperfine structure (due to the nuclear magnetic moment); 

• recoil corrections; 

• radiative (QED) corrections; 

• nuclear- structure corrections. 

A structure of levels with the same value of the principal quantum number n is 
a kind of signature of any atomic system. For most of the precision applications 
the substructure of interest is related to n — 2. The corrections decrease with 
a value of the principal quantum number as n~ 3 or faster. The only exception 
is the Uehling correction for muonic and exotic atoms which scales as rC 2 for 
medium and high Z . 
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Table 1 

Various contributions to the energy levels. The results are in units of (Za) 2 mc 2 , 
where m is the mass of the orbiting particle. Here: M is the nuclear mass and 
m p is the proton mass which enters equations if one measure the nuclear magnetic 
moment in units of the nuclear magneton. A contribution of the nuclear magnetic 
moment, i.e., the hyperfine structure, appears if the nuclear spin is not zero. i?jv 
stands for the nuclear (charge) radius. 

We summarize the most important corrections to the levels determined by 
Eq. (3) in Table 1. The main phenomena essentially contributing to the struc- 
ture of the energy levels in hydrogen-like atoms are listed below. 

• Non-relativistic consideration 

- The dominant contribution is determined by the Schrodinger equation for 
a particle with mass m bound by a central Coulomb field. The result has 
order (Za) 2 mc 2 and we discuss all other corrections in units of (Za) 2 mc 2 . 

- The leading non-relativistic correction is already incorporated into Eq. (3) 
by the introduction of the reduced mass (4). The correction has a frac- 
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tional order of m/M, and the result, presented in Eq. (3), is already cor- 
rected. This correction is responsible, for example, for the isotopic shift of 
hydrogen and deuterium gross structure. In electronic atoms, the correc- 
tion has order m e /M ~ 5 ■ 10~ 4 /A, where A is the atomic number (i.e., 
the atomic mass in the universal atomic mass units). The isotopic shift 
between atoms with atomic numbers A and A' = A + AA involves for 
heavy elements an additional small factor of AA/A. In muonic atoms the 
correction is of order of 0.1/ A and in positronium the effect reduces the 
result in the leading order by a factor of two. 
The result (3) depends on a principal quantum number n only. The structure 
of levels due to the non-relativistic approximation does not have a single 
common notation. It used to be named as the gross structure, the main 
structure, the Rydberg structure, the basic structure etc. 
• Relativistic effects 
- The leading relativistic corrections are determined by a solution of the 
Dirac equation 

E D = mc 2 x F D (nlj) , (5) 



where 



F D (nij) = (l + 



(Zaf 



(n-j -1/2+ J(j + 1/2)2 _ (Zq)2 ) ' 



-1/2 



(6) 



and thus 



E n = mc 2 x 



_ {Zaf _ {Zaf_ ( 1 _ 

2n 2 2n 3 I j + 1/2 An) 



(7) 



The first term above is for the rest energy, the next is the leading non- 
relativistic term and the third term (of order of (Za) 4 mc 2 in absolute 
units, or [Za) 2 in units of (Za) 2 mc 2 which we apply in this section) is the 
relativistic correction. This relativistic correction splits levels with the 
same n, but different values of the angular momentum j. However, the 
result does not depend on the orbital momentum /. In particular, these 
relativistic effects are responsible for the fine structure in the hydrogen 
atom. They split the 2p 3 / 2 and 2p 1 / 2 states, but the energy of the 2si/ 2 
and 2p 1 / 2 is still the same. 

The Dirac equation itself is an equation ignoring the nuclear motion. Both 
the nonrelativistic reduced-mass effects and the leading relativistic recoil 
correction, which contributed in order (Za) 2 m/M, can be expressed in 
terms of the Dirac energy (see, e.g., [15]) 



AE TOC = mc 2 x 



rriR — m 
m 



(F D -l)- 



m 



R 



2Mm 2 V 



(F D -1 



(8) 
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We note that a simple substitution m — > ttir reproduces the correct results 
for the leading terms for both gross and fine structure since only the first 
term in (8) contributes there. However, the substitution leads to an incorrect 
result for the leading relativistic correction for, e.g., the Is state, where the 
second term in (8) is also important. 

• Contributions involving the nuclear structure effects appear already in 
the non-relativistic approximation. 

- One of the main effects is the hyperfine structure, a result of interaction 
of the nuclear spin / with the angular momentum of the orbiting particle 
j. The energy of the levels depends on the quantum number F (where 
F = j+I). The result is of order of (Za) 2 m/M if we use the notation of the 
particle physics, where considers a nucleus as a particle of a charge Ze and 
mass M and measures its magnetic moment in units of its own magneton 
Zeh/2M. Here e is the proton charge. This is the case of positronium and 
muonium. 

In the case of a compound nucleus (deuterium, tritium, helium-3 etc.) 
another convention is in use. The magnetic moment is measured in units of 
the nuclear magneton eh/2m p , where m p is the proton mass. In such a case 
one usually speaks about a hyperfine structure of order of a(Za)m/m p . 

Indeed, these two notations are related to the same order of magnitude 
since Z/M ~ Z/(Am p ) ~ l/m p with Z < A < 3Z. In the case of the 
hydrogen atom these two notations are completely identical. 

The hyperfine structure has the same parametrical order for states with 
any I and j, however, the result for the s states is usually considerably 
larger numerically. 

A value of the nuclear magnetic moment is a consequence of its nuclear 
structure and in conventional atoms (hydrogen, deuterium etc.) cannot be 
calculated ab initio in contrast to muonium and positronium. 

- In the leading order, effects of a distribution of the nuclear charge (so- 
called finite-nuclear- size effects) shift energy levels of the s states only 
and the correction has order {ZamcRjss /h) 2 , where is a characteristic 
nuclear radius. 

• Quantum Electrodynamics offers a variety of effects and their hierarchy 
depends on a kind of atom. 

- The leading QED effect in the hydrogen atom (and other conventional 
atoms) is the self energy contribution (see Fig. 2) which leads to a correc- 
tion of order a(Za) 2 ln(l/(Za)). It splits the levels with the same j, but 
different values of I. An example is a splitting of 2si/ 2 — 2pi/ 2 (the Lamb 
shift). 

- The vertex diagram in presence of the magnetic field is chiefly responsible 
for a correction to the magnetic moment (i.e., for the anomalous magnetic 
moment of an electron) and other effects are of higher order. As we al- 
ready mentioned, the anomalous magnetic moment of an electron was first 
detected as an anomaly in the hyperfine structure of atomic hydrogen. 

- An electron vacuum polarization (see Fig. 3) in the hydrogen atom is 
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Fig. 2. The self-energy diagram which is responsible for the dominant contribution 
to the splitting of 2s]y 2 — 2pi/2 (the Lamb shift) in the hydrogen atom. The doubled 
electron line is related to the Green function of an electron in the external Coulomb 
field. 

responsible for a small fraction of the Lamb shift. In a general case, the 
correction, also called the Uehling correction, has order aF(Zam/m e ), 
where 







x < 1 


F(x) ~ < 


1, 


x ~ 1 




In a; , 


X > 1 



In the conventional atoms, where m e = m and x = Za, the Uehling 
contribution does not change the hierarchy of the intervals. In muonic 
atoms, where x ~ 1.5 Z, the Uehling correction can produce 2s — 2p 
splittings bigger than 2p 3 / 2 — 2pi/2- The hierarchy depends on Z. 

I 

* 

Fig. 3. The vacuum polarization diagram (the Uehling contribution). In the light 
muonic atoms it leads to the dominant contribution to the splitting of 2s — 2p (the 
Lamb shift) at low and medium Z. 

- In the case of positronium the two effects above are not so important as 
the virtual annihilation of the bound electron and positron into a single 
photon (see Fig. 4). The one-photon annihilation splits two hyperfine levels 
related to the S states: the 3 Si state (triplet) is shifted while the 1 S state 
(singlet) is not affected. 

- The QED effects are also responsible for decay of atomic states (see Fig. 5). 
The radiative decay line width is basically a(Za) 2 . This result is related to 
the one-photon electric-dipole transition (so called El transition) which is 
a dominant decay mode for all levels in light hydrogen-like atoms except 
the 2s and Is states. The upper hyperfine component of the Is states can 
decay via a magnetic-dipole (Ml) transition with an extremely small line 
width while the metastable 2s state decays via a two-photon electric-dipole 
(2E1) transition with the line width of order of a 2 (Za) 4 . 
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Fig. 4. The one-photon annihilation diagram. The real annihilation into a single 
photon is not possible for kinematic reasons, while the virtual one is possible and 
it shifts the energy levels of s states with total spin 1. 



Fig. 5. The imaginary part of the self energy diagram is the energy width of the 
level and thus it is equal to the probability of its radiative decay via emission of a 
single photon. 

- In the case of positronium the line width of the S states 1 (n = 1,2) 
depends on their total angular momentum (spin). The triplet states ( 3 5i) 
are long living decaying via the three-photon annihilation, while the two- 
photon annihilation is a dominant mode for the short-living singlet states 
( 1 S'o). The widths are a 4 and a 3 respectively. Higher S states additionally 
to annihilation modes can also decay radiatively into P states (one photon 
El transition) with a related width of order of a 3 . In the case of non-S* 
states (P, D etc.) the radiative decay dominates. 

In any case the contribution in Eq. (3) dominates for an interval between 
levels with An ^ 0. However, the structure of levels with the same value of 
n depends on details of the atom. In Fig. 6 we present three different basic 
spectra of the structure at n — 2. 

• The first one is realized in 'normal' (electronic) hydrogen-like atoms (hy- 
drogen, deuterium, helium ions etc.). The muonium spectrum is of the same 
kind. The largest splitting, of order (Za) 4 m e c 2 , is the fine structure (i.e., a 
splitting between levels with a different value of the electron angular mo- 
mentum j), the Lamb shift arising from the electron self-energy effects is of 



1 We use the low case letters (s, p etc.) for single-electron atoms, while the capital 
letters are used for atoms with a few electrons. Other indexes for single-electron 
atoms are related to the electron properties, while a total angular momentum F is 
indicated in the parentheses: e.g., ls 1 / 2 (-P = 1). In the case of positronium, since the 
masses of electron and nucleus are the same, we use capital letters indicating that 
we deal with the total atomic momentum related to two particles and the subscripts 
and superscripts are related to properties of an atom as a whole. 
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4 He), while others have a non-zero spin (/ = 1/2 in hydrogen, 
helium-3 and I = 1 in deuterium etc.). In the latter case, the 



Fig. 6. Scheme of the lowest excited levels (n = 2) in different simple atoms (not to 
scale) . 

order a(Za) 4 m e c 2 In (l/(Za)) and it splits the levels with the same j and 
different values of the electron orbital momentum I. Some nuclei are spinless 
(like, e.g. 
muonium 

interaction with the nuclear spin splits levels with the same electronic quan- 
tum number. The splitting is of order [Zot) 4 w? e c 2 /M or a(Zaf'm 2 c 2 /m p , 
where M is the nuclear mass, and the structure depends on the value of the 
nuclear spin. The scheme in Fig. 6 is for a nuclear spin of 1/2. A difference 
for higher I is that the states with j > 1 (such as 2p 3 / 2 with j = 3/2) 
are split in three or more subcomponents. In a general case, the number of 
hyperfine components is the smaller number among 2j + 1 and 21 + 1. 

• The structure of levels in positronium and muonic atoms is different because 
other QED effects enter consideration. For positronium, an important fea- 
ture is a real (into two and three photons) and virtual (into one photon) 
annihilation. The former is responsible for the decay of the s-states, while 
the latter shifts triplet s-levels (l 3 Si and 2 3 Si in particular). The shift is of 
the order of a 4 m e c 2 . Contributions of the same order arise from relativistic 
effects for the both, the electron and the nucleus (namely, positron), and 
from hyperfine interactions. As a result, the structure of the positronium 
levels at n = 2 has no apparent hierarchy (Fig. 6). 

• Another situation is for the spectrum of the muonic atoms. A difference to 
the hydrogen case comes from a contribution due to the vacuum polarization 
effect (the Uehling potential). Effects of electronic vacuum polarization shift 
all levels to the order of a(Z«) 2 m M c 2 . This shift is a nonrelativistic one and 
it splits 2s and 2p levels, but in does not depend on j and cannot split 
2^3/2 and 2pi/2- The fine and hyperfine structures are of the same form 
as for the normal atoms (i.e., {Za) 4 m^c 2 and (Za) 4 m^c 2 /M, respectively), 
but the latter is much more important than in a conventional atom since 
m fj,/ m p ~ 1/9. At the case of low Z the Lamb shift induced by the Uehling 
potential is the dominant correction to the energy levels. 

Table 1 shows Z dependence of various crucial contributions to energy of 
bound states in two-body atoms and to their decay width. Since the radiative 
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linewidth increases with Z, it is much easier to work with light atoms which 
provide a longer lifetime of excited states. We note, however, that importance 
of QED effects also increases with Z in such a way that a ratio of a correction 
to width, which is of interest, rather favors heavy ions with high Z. Meanwhile 
there is a number of other reasons which make \ow-Z atoms more attractive. 

• Studies of neutral atoms and ions involve quite different methods and as a 
result, some quantities are better determined for neutral two-body atoms 
(like, e.g., the hyperfme interval in the ground state) while others lead to a 
higher accuracy for ions (like, e.g., the hyperfme interval for the excited 2s 
state in helium-ion 3 He + ). Indeed, all neutral two-body atoms have a low 
value of the nuclear charge, namely, Z — 1. 

• Precision QED tests are realized in atomic systems which can be calcu- 
lated ab initio. Such a theory can be in principle developed for any QED 
effects, but not for the nuclear structure. Importance of the nuclear struc- 
ture effects increases with Z. In the leading non-relativistic approximation 
(see Table 1) the related nuclear-spin-independent contribution to the en- 
ergy is proportional to Z 4 and a squared value of the nuclear charge radius 
which also increases with Z. In the case of the nuclear-spin-dependent con- 
tributions (hyperfme effects) the related correction is proportional in the 
non-relativistic approximation to Z 5 and to a characteristic value of the nu- 
clear size. In the case of high Z the non-relativistic approximation still leads 
to a correct order of magnitude, however, the dependence on the nuclear 
charge and nuclear radii becomes more complicated. 

As a consequence, studying hyperfme interactions one has either to study 
the nucleon-free atoms (muonium, positronium) or try to cancel the nuclear 
effects between a few measured quantities. That is possible because low- 
Z atoms allow to present the nuclear effects in terms of a few effective 
parameters (such as the nuclear charge radius, nuclear magnetic radius etc.). 
On contrary, for highly charged ions we rather need to apply certain models 
of distribution of their charge and magnetic moment and thus the QED 
tests with high-Z ions become model-dependent, in comparison with model- 
independent calculations for low-Z. 

• The Rydberg constant is related to an ultraviolet part of spectrum. Thus, 
at low Z (Z — 1,2) there is a good chance to study some transitions in gross 
structure with lasers while for higher Z laser spectroscopy can be used only 
for the fine structure and the Lamb shift as well as for muonic atoms. Strong 
dependence of the transition frequency on Z allows only a few opportunities 
for a spectrum of two-body atoms with Z ^ 1 and/or m ^ m e . 

• In contrast, a number of transitions in neutral atoms Z — 1 are suited 
for laser spectroscopy. That is not a surprise. The characteristic atomic 
and molecular energies we meet in our life are of the same order as in 
hydrogen since they are mostly related to certain single-electron excitations 
in a Coulomb field of an effective charge of Z e s ~ 1. Since we are forced 
to meet them in our common-day, industrial and scientific life the most 
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advanced radiation sources and useful tools were developed for these regions 
(optical with some inclusion of infrared and ultraviolet domains). Their 
development has a history of a few centuries of making lenses, mirrors, 
prisms etc. 

• As we noted, the ratio of the crucial QED corrections of interest (e.g., 
a 2 (Za) 6 mc 2 contributions to the Lamb shift) and the lifetime (which scales 
for the El transitions such as 2p — > Is as a 5 (Za) A mc 2 ) clearly favors higher 
values of Z. However, that is not a whole story. First, there is a metastable 
2s state in hydrogen-like atoms which lifetime scales as a 2 (Za) e mc 2 and use 
of this level is crucial for a number of experiments. Because of the longer 
lifetime and a proper value of the Is — 2s interval it is much easier to deal 
with this level at low Z. In the hydrogen atom the 2p lifetime is about 
10% of the Lamb splitting 2s±/ 2 — 2p ± / 2 while a much smaller value of the 
natural 2s lifetime cannot be even seen from the line shape of the two- 
photon Is — 2s excitation in precision spectroscopy experiments because of 
various perturbations which significantly shorten the 2s lifetime (which is 
still much longer than that of the 2p state). The narrow Is — 2s transition 
was successfully used in a number of hydrogen experiments while there is no 
access to this narrow transition for high Z. Secondly, the lifetime of the El 
width sharply depends not only on a value of the nuclear charge Z but also 
on the principal quantum number n. With such an n-dependence which can 
be roughly approximated as 1/n 3 one can realize that the 2p transition is 
one of the broadest. The high-Z experiments used to deal with 2s — 2p and 
Is — 2p intervals and thus the involvement of the 2p state is unavoidable. 
On contrary, the recent success in hydrogen spectroscopy is related to high 
n levels studying 2s — nd and 2s — ns at n = 8 — 12. Use of the metastable 
2s level and highly excited levels in hydrogen compensates all disadvantages 
of the correction-to-width ratio for low Z and makes the low-Z case most 
favorable for the bound state QED tests. 

Since experimental and theoretical methods favor the light atoms, we focus 
our consideration of particular QED tests on hydrogen and other low Z atoms 
and touch other atoms only in case if they provide any results competitive 
with Z = 1 atomic systems. 



3 Optical measurements in hydrogen atom and determination of 
the Rydberg constant 

About fifty years ago it was discovered that in contrast to the spectrum pre- 
dicted by the Dirac equation, there are some effects in hydrogen atom which 
split the 2si/2 and 2pi/ 2 levels. Their splitting known now as the Lamb shift 
(see Fig. 7) was successfully explained by quantum electrodynamics. The QED 
effects lead to a tiny shift of energy levels and for thirty years this shift was 
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studied by means of microwave spectroscopy (see, e.g., [31,32]) measuring ei- 
ther directly the 2sx/ 2 — 2py 2 splitting or the larger 2p 3 / 2 — 2si/ 2 splitting (fine 
structure) where the QED effects are responsible for approximately 10% of 
the interval. 

^^^^^^^^Hl 2 P3/2 

rf transition (fine structure) 

2s-|/2 a ~ rf transition (Lamb splitting) 

2P1/2 



1s 1/2 . 



two-photon uv transition (gross structure) 



Fig. 7. An example of different transitions in the hydrogen atom (not to scale). 
Transitions within the fine structure and Lamb splitting are in the microwave range, 
while the Is — 2s two-photon transition lies in the ultraviolet domain. The Lamb 
splitting is the difference of the Lamb shifts of 2s\/2 and 2pi/2 levels. The hyperfine 
structure is neglected. The 2p linewidth is approximately 10% of the Lamb shift, 
which itself is approximately a tenth part of the fine structure splitting. 

A recent success of two-photon Doppler-free spectroscopy [33] opens another 
way to study QED effects directed to high-resolution spectroscopy of gross- 
structure transitions. The energy of such a transition between levels with 
different values of the principal quantum number n is mainly determined by 
the Coulomb-Schrodinger formula in Eq. (3). However, an unprecedently high 
accuracy achieved by to date allows to study various tiny perturbations to 
Eq. (3), caused by effects of relativistic physics, quantum electrodynamics 
and nuclear structure. All studied transitions are related to the ultraviolet 
part of the spectrum. 

For any interpretation of hydrogenic lines in terms of QED effects one has to 
determine a value of the Rydberg constant 

„ a 2 mc , . 

*~ = ^r- < 10 > 

One more problem for the interpretation of the optical measurements is the 
involvement of few levels, significantly affected by the QED effects. In con- 
trast to radiofrequency measurements, where the 2s — 2p splitting was under 
study, the optical measurements have been performed with several transitions 
involving a number of states (Is, 2s, 3s etc.). It has to be noted that the Lamb 
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shift for levels with I ^ has a relatively simple theory, while theoretical cal- 
culations for the s states lead to several serious complications. The problem 
has been solved by introducing an auxiliary specific difference [34,35] 



A(n) = E L (ls) -n 3 E L (ns) 



for which theory is significantly simpler and more clear than for each of the s 
states separately 

Combining theoretical results for the difference [35,36] with a measured fre- 
quency of two or more transitions one can extract a value of the Rydberg 
constant and the Lamb shift in the hydrogen atom. 
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Fig. 8. Progress in determination of the Rydberg constant by means of two-photon 
Doppler-free spectroscopy of hydrogen and deuterium. The label CODATA stands 
for the recommended value of the Rydberg constant i?oo(1998) [21] from Eq. (12). 
The most recent original value is a preliminary result from MIT obtained by mi- 
crowave means [37]. 

The recent progress in determination of the Rydberg constant is presented in 
Fig. 8 (see Refs. [33,21] for references). The progress is also clearly seen from 
the recommended CODATA values of 1986 [20] and 1998 [21]: 

#00(1986) = 10 973 731.534(13) m' 1 



and 



^(1998) = 10 973 731.568 549(83) m" 1 



(12) 



The former value was derived from one-photon transitions (Balmer series) and 
was 'slightly' improved later (by a factor of 4.5), but all further progress that 
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led to the 1998's value (improvement by more than 30 times) was a result of the 
study of two-photon transitions in hydrogen and deuterium. Figure 8 shows a 
comparison of several recently published values for the Rydberg constant. 

As mentioned before, the other result, which can be extracted from the optical 
spectroscopy of the hydrogen atom, is the Is (2s) Lamb shift. However, studies 
of the Rydberg constant and the Lamb shift, which are related quantities, 
are still not the same. While the optical measurements deliver us the most 
accurate data on both values (and in fact the data are strongly correlated), 
the microwave experiments still supply us with various less precise data on 
both quantities and the microwave data on the Lamb shift and the Rydberg 
constant are completely independent from each other. 

The only example of a microwave determination of the Rydberg constant is 
an MIT experiment [37], which dealt with highly excited hydrogen levels. The 
transitions under study were between n = 30 and n = 27 for high / and because 
of such high values of n and / the levels were not very sensitive to the QED 
effects and thus the derived value of the Rydberg constant is 'Lamb-shift-free'. 
The preliminary result [37] 

#00 = 10 973 731.568 34(23) m" 1 

agrees with other values (see Fig. 8). All other important related microwave 
results in the field have been achieved for either the Lamb shift or the fine 
structure. 



4 The Lamb shift in the hydrogen atom 

After the progress of the last decade, the optical determination of the Lamb 
shift [33] superseded traditional microwave measurements [31,32], but still 
there is a number of compatible microwave results among the published data. 
The experimental results are summarized in Fig. 9. 

To discuss the state of the art in study of the Lamb shift in detail let us start 
with the experimental data. A number of different intervals in the hydrogen 
atom have been precisely studied for about a century. Bound state QED and in 
part even QED itself were established essentially in order to explain the Lamb 
shift and the fine and hyperfine structure in the hydrogen atom as well as in the 
deuterium atom and the helium ion 2 . The corresponding transitions lie in the 

2 The Lamb shift is sometimes named after Lamb and Retherford, referring to their 
famous papers on the Lamb shift in the hydrogen atom [40]. However, that is not 
quite correct, since at approximately the same time a result on the Lamb shift in 
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Fig. 9. Measurement of the Lamb shift in the hydrogen atom. The most accurate 
experimental result comes from a comparison of the Is — 2s interval measured 
at MPQ (Garching) [38] and the 2s — ns/d intervals at LKB (Paris) [39], where 
n = 8, 10, 12 (see also [33] for detail). Three more results are shown for the average 
values extracted from direct Lamb shift measurements, measurements of the fine 
structure and a comparison of two optical transitions within a single experiment 
(i.e., a relative optical measurement). The filled part is for theory. Theory and 
evaluation of the experimental data are presented according to Ref. [36]. 

microwave part of the spectrum. In the last decades, progress in radiofrequency 
experiments with hydrogen and other light atoms was rather slow. The results 
of over than thirty years related to the Lamb shift in hydrogen are presented 
in Fig. 10 (measured directly) and Fig. 11 (the Lamb shift deduced from the 
measured fine structure interval 2p 3 / 2 — 2s! / 2 ) . The most recent results are 
presented separately, while the older experiments are averaged. The highest 
accuracy was claimed in [45] and we correct here their value according to 
[34,46]. Nevertheless, the result is not included in the Lamb shift average 
value in Fig. 9 (the LS value) because of unclear status of its uncertainty, 
which is discussed in part below in this section. 



To reach a value of the Lamb shift from the fine structure (2p 3 / 2 — 2si/ 2 ) 
measurement we need to use a value of the 2p 3 / 2 — 2pi/ 2 splitting which was 
found theoretically 



AE d 2 Pj ) = AE L (2p) + x [a c ■ A %> (^ + C R 



a(Za) 2 



An In + + G SE (Za 



7r y (Za) 
+ higher order terms , (13) 



the helium ion was presented in a paper by Skinner and Lamb [41]. 
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Fig. 10. Direct measurement of the Lamb shift in the hydrogen atom. The most 
recent result was obtained in [42,43], while the older results [44] are averaged. The 
result of Sokolov and Yakovlev [45] has been corrected according to [46]. 



Hagley & Pipkin 






i ♦ 1 


Older 


i 


♦ 


measurement 



1 057 800 kHz 1 057 880 kHz 



Fig. 11. An indirect determination of the Lamb shift (2s]/2 — ^Pi/2) m atomic 
hydrogen via a study of the fine structure 2p 3 / 2 ~ ^ s i/2- The most recent result was 
obtained in [47], the older results [48] are averaged. 

where we introduce a j-independent part of the Lamb shift of the 2p state 



AE L( 2 P ) = A_L m J (_5) |--ln fco (2p)x (l + — ) 

+ -V- 111 2 + ir ( c - + 3 {Za) ) } 

+ higher order terms . (14) 

The state-dependent coefficients A2 , v4 61 , Aq^, CR ec , Cr and B 62 and the 
function Gse{Zo) are well known (see in Table 2); a e stands for the anomalous 
magnetic moment of the electron and lnfco(2p) = —0.030 016 709... is the 
Bethe logarithm. 

The direct measurements of the Lamb splitting (i.e., the interval between 2s\/2 
and 2^1/2) need in principle no auxiliary QED calculations. An exception is 
the result of experiment [45] where the originally measured quantity was a 
product of the 2si/ 2 — 2p 1 / 2 splitting and the lifetime of the 2p 1 / 2 state. The 
measurement was claimed to be the most accurate determination of the Lamb 
shift (the assigned uncertainty was about 2 ppm). The result [45] has been cor- 
rected due to a recalculation of the lifetime [34,46], where the leading radiative 
correction was taken into account within a logarithmic approximation 
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Table 2 

Coefficients for higher-order contributions to the hydrogen Lamb shift AEi(2p) and 
specific difference A(2) (see Eq. (11)). Most of the coefficients are universal and are 
the same for, e.g., deuterium, with one exception: coefficient Cr depends on the 
nuclear spin, it is presented for hydrogen (/ = 1/2); coefficient C m j x = (g p ) 2 /162 
is due to mixture of the 2pi/ 2 (F = 1) and 2p^/ 2 {F = 1) states by the hyperfine 
interaction and was considered in [49] (see, also [50]) and F is the complete angular 
momentum and its value for the 2p 1 / 2 state in hydrogen can be either or 1, while 
for the 2p 3 / 2 state it can take value of 1 or 2. 
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We remind here that Roo is related to the wave number iy/c), while oRoo 
corresponds to frequency v. Factor of 2rr appears because hu = Hcu and the 
lifetime r is related to the line width of the angular frequency cu = 2itv. 

The logarithmic part of a(Za) 2 term was calculated in [34,46] (see also [51,52]). 
Recently a non-logarithmic term was also calculated C = 6.576 . . . [53]. Since 
our result [52] already included an estimation of the constant Co = 6.2(5), that 
leads to a marginal shift in the value of Co, but makes the result previously 
based on a naive estimation [52] more reliable. We also have to mention that 
an exact calculation of C alone is not sufficient for any application because 
the line width is not well defined. As it is well known, the line shape contri- 
butions to the width enter in fractional order h/[r ■ AE], which in the case of 
an El transition in a hydrogen-like atom is related to an a(Za) 2 correction. 
That is of the same order as under question and thus the non- logarithmic term 
should depend on the phenomenon under study. If we like to express the result 
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of the Sokolov's experiment [45] in terms of T{2p 1 / 2 ) ■ [El(2si/ 2 ) — El(2pi/ 2 )} 
and to present the effective lifetime r(2pi/ 2 ) in terms of Eq. (15), we see that 
the coefficient Cq is phenomenon-dependent. However, we note that a numer- 
ically big contribution to the effective constant C cannot appear in such a 
way since %/[r{2p) ■ AE(2p — Is)] ~ 0.33 a 3 / 71 " • To be conservative we double 
the value of h/[r(2p) ■ AEi{2p — Is)] and estimate the phenomenon-dependent 
contribution to T(2p 1 / 2 ) as ACo — ±0.7 which is to be added to Co for the 
interpretation of experiment [45]. 

We have not taken into account the result of [45] when calculating the average 
values in the summary picture (Fig. 9) because there are certain doubts in 
its accuracy noted by Hinds [31] and it is not really clear if the result is as 
accurate as it was stated originally. The result of [45] is rather in a not good 
agreement with both the theoretical value and a grand average over the other 
data. The discrepancy with theory does not exceed two standard deviations 
and it is too early to consider that as a real contradiction. 

We need, however, to note that a common opinion on the future of the mi- 
crowave measurements of the Lamb shift in hydrogen and on the reliability 
of the result [45] with a 2-ppm uncertainty has assumed two contradicting 
statements. 

• First, it is generally believed that since any 2s — 2p intervals, (either 2s\j 2 — 
2pi/2 (the Lamb splitting) or 2p%/ 2 — 2si/ 2 (the fine structure)) have a width 
of 0.1 GHz (because of the radiative lifetime of the 2p state), the Lamb 
splitting (which is about 1 GHz) cannot be measured better than at the 
10 ppm level. Thus, it is expected that the statistical error of a microwave 
measurement should be larger than 10 ppm. 

• Secondly, sharing the Hinds' point of view [31], it is believed that the accu- 
racy of the Sokolov's experiment [45] is not as high as claimed only because 
of a possible systematic error in the determination of the atomic beam ve- 
locity, which, in principle, might be somehow fixed in the next generation 
of experiments. Meanwhile, there has been no doubt expressed on the sta- 
tistical treatment of the Sokolov's data. 

Because of this inconsistency in a opinion of experts, the question on a possi- 
bility of a further progress of microwave measurements of the Lamb shift with 
an uncertainty below 10 ppm remains unclear. A measurement of the Lamb 
shift in deuterium within the Sokolov's scheme could provide a chance to test 
partly systematic effects since the result may be compared with the hydrogen- 
deuterium isotopic shift of the Is — 2s interval [54]. The latter, combined with 
the theoretical result on specific difference in Eq. (11) [34,35,36], leads to the 
isotopic shift for the Lamb shift 
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with an uncertainty below 1 kHz. Here a calculation of three terms (namely, 
the Lamb shift of the 2pi/ 2 state, the specific difference of Eq. (11) at n = 2 
and the energy related to an effective Dirac equation 



2Mm v 

•_>//•• M 2 m 2 U + l/2 / • 1 2 1 " 



(17) 



which takes into account recoil effects up to (Za) 4 mlc 2 /M 2 (the dimension- 
less Dirac energy Fo(nlj) is defined in Eq. (6)), can be performed with an 
uncertainty substantially below 1 kHz. The only experimental value in the 
right-hand part of the equation (the isotopic shift of the Is — 2s frequency) is 
also known with an uncertainty below 1 kHz [54]. 

Essential progress in studies of the hydrogen Lamb shift came some time ago 
from the optical two-photon Doppler-free experiments (see [33,55] for detail). 
The Doppler-free measurements offer a determination of some transition fre- 
quency in the gross structure with an accuracy high enough to use the results 
to find the Lamb shift. However, the cancellation of the linear Doppler effect 
is not an only essential advantage of those experiments. A microwave mea- 
surement is to be designed to determine an interval between the 2s state and 
either the 2pi/ 2 or 2p 3 / 2 state which has a broad line width limiting the accu- 
racy of the measurement. It is not absolutely clear what the ultimate limit of 
the measurements involving the 2p states is, but indeed an experiment with 
more narrow levels should have a good chance for a higher accuracy. That is 
the two-photon excitation of the metastable 2s state to higher excited ns and 
nd levels that allows to deal with some relatively narrow levels because of the 
n~ 3 dependence of their line width. However, with the expanding number of 
the involved energy levels, two special problems as we mentioned before have 
appeared in the optical experiments: a determination of the Rydberg con- 
stant and needs to reduce the variety of the QED contributions for a number 
of levels to a single quantity, e.g., to the Lamb shift of the 2s (or Is) state 
E L (2s). 

Two methods were applied to manage the problem with the Rydberg contri- 
bution. The first of them was for a measurement of two different frequencies 
within one experiment with the ratio of the frequencies being approximately 
an integer number. Extracting a beat frequency one can avoid the problem of 
determining the Rydberg constant at all. Three experiments have been per- 
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formed in this way: the Garching experiment dealt with the Is — 2s transition 
and the 2s — 4s (and 2s — Ad) transition [56], at Yale the Is — 2s frequency was 
compared with the one-photon 2s — 4p transition [57] and actually that was 
the only precision optical experiment with an one-photon transition among 
the recent generation of the experiments on hydrogen spectroscopy. More re- 
cent Paris experiment dealt with the Is — 3s and 2s — 6s (and Qd) intervals 
[58]. The values derived from those experiments are collected in Fig. 12. 



Comparison of 1s- 2s (Garching) with 2s-8s/d (Paris) 

I 8 1 



ROM (Paris) 
• ROM (Garching) 

• ROM (Yale) 

1 057 830 kHz 1 057 870 kHz 

Fig. 12. An optical determination of the Lamb shift in the hydrogen atom. The label 
ROM stands for relative optical measurements, which dealt with a comparison of 
two optical frequencies whithin the same experiment. 

Another way to manage the problem with a determination of the Rydberg 
constant is to perform two independent absolute frequency measurements (i.e., 
the measurements with respect to a primary cesium standard). Comparing 
them afterwards one can determine both the Rydberg constant and the Lamb 
shift. Combining results from Garching on the Is — 2s transitions [38] and from 
Paris on the 2s — 8s/8d/12d intervals [39], the most accurate optical value of 
the Lamb shift has been reached (Fig. 12). 

Some of the optical experiments were also performed for deuterium (see, e.g., 
[39]) and that may improve the accuracy in the determination of the Rydberg 
constant and, thence, indirectly, of the hydrogen Lamb shift. 

The recent progress with the visible and ultraviolet optical measurements 
became in part possible because of use of the frequency comb generator, an 
advanced technique which allows to compare two optical frequencies to each 
other and an optical frequency to a microwave frequency [59]. Because of 
the frequency comb generation, a comparison of two transitions within one 
experiment gives no advantages anymore. We expect that further progress 
will be due to independent absolute measurements of different transitions. 

The optical experiments delivered data which involve a number of levels (Is, 
2s, 3s etc.). Therefore, a new problem appeared in the increasing number of 
levels and each of these levels has its own Lamb shift, the value of which should 
be determined somehow. A large number of quantities to extract from exper- 
imental data may reduce advantages of the optical approach. A problem with 
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the number of Lamb shifts was solved with the help of the specific difference 
in Eq. (11) which is found to be of the form [34,35] 
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where the n-dependent coefficients A 6 i, A^, C Rec and B 62 are presented for 
n = 2 in Table 2. The coefficients for other n and a table for the function 
G (Za) can be found in Refs. [35,36]. The uncertainty was also discussed 
there. A crucial issue is that the difference has a better established status than 
that for the Is (or 2s) Lamb shift. In particular, various QED corrections van- 
ished for the difference, some other are of a reduced numerical value. The 
finite-nuclear-size contribution, which produces the largest theoretical uncer- 
tainty, is cancelled out after combining the Is and 2s contributions. 

To analyze the experimental data, we need to use a certain piece of the bound 
state QED theory, such as results on the Lamb shift of 2pj states given in 
Eqs. (13) and (14), on the specific difference A(n) (18) defined in Eq. (11) 
and a result for the lifetime of the 2p state (15). 

Thus, prior to test the QED theory of the hydrogen Lamb shift, some the- 
oretical results have to be applied to extract the experimental values of the 
Lamb shift from the data of the measurements. Due to that we have to clarify 
here the word 'theoretical'. To our mind, a value is a theoretical one if it is 
sensitive to theoretical problems (like, e.g., the problem of the determination 
of the fundamental constants and phenomenological parameters needed as an 
input to the calculations and the problem of a proper estimation of uncalcu- 
lated higher-order QED corrections for the Lamb shift). A calculated value, 
which is not sensitive to this kind of problems, is not theoretical, but rather 
a mathematical one. 

Rigorously speaking, there is neither pure theoretical value nor experimental 
value. Theoretical calculations involve experimentally determined values of 
certain parameters, or even functions, while any measurement needs a theory 
to interpret measured quantities in terms of collected records of experimental 
devices. 

The main problem of theoretical calculations for the hydrogen Lamb shift is 
a lack of an accurate determination of the proton charge radius and an in- 
sufficient understanding of uncertainty due to the higher-order QED effects. 
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The former source of theoretical uncertainty, related to experimental particle 
physics, dominates. The proton charge distribution does not affect evalua- 
tion of experimental data via E L (2pj), A(n) and r(2p), while the theoretical 
uncertainty due to the higher-order QED effects enters the evaluation. 

An additional discussion to clarify the problem is needed, since the auxiliary 
piece of theory needed to evaluate the experimental data for the Is Lamb shift 
is related to the Lamb shift theory, which in principle we would like to verify 
measuring the Is Lamb shift. A question on self-consistency of such a test 
of the QED theory of the Lamb shift arises. We emphasize, however, that a 
crucial point is that the specific difference A(n) for the ns states, which are 
involved in actual experiments, and the Lamb shift and the lifetime of the 
2pj states corresponds to a part of theory which is understood substantially 
better than theory of the Is Lamb shift. The higher-order contribution enters 
auxiliary theoretical expressions for E L (2pj), A(n) and r{2p) in an essentially 
weakened form, since the dominant contributions are proportional to <5/o- 

A successful deduction of the Lamb shift (Fig. 9) in the hydrogen atom pro- 
vides us with a precision test of bound state QED and offers an opportunity 
to learn more about the proton size. Bound state QED is quite different from 
QED for free particles. The bound state problem is complicated by itself even 
in the case of classical mechanics. The hydrogen atom is the simplest atomic 
system; however, a theoretical result for the energy levels is expressed as a 
complicated function (often a perturbation expansion) of a number of small 
parameters [60] (see review [16] for a collection of theoretical contributions): 

• the fine structure constant a, which counts the QED loops; 

• the Coulomb strength Za\ 

• the electron-to-proton mass ratio; 

• the ratio of the proton size to the Bohr radius. 

Indeed, in the hydrogen atom Z — 1, but it is customary to keep Z even 
for hydrogen in order to trace out the origin of the corrections, especially 
because the behavior of expansions in a and Za differs from each other. In 
particular, the latter involves large logarithms (\n(l/Za) ~ 5, \n 2 (l/Za) ~ 24 
and ln 3 (l/Zo;) ~ 120 at Z = 1) [60,61] and big coefficients. Comparing a and 
Za we note that in the relativistic case, when the energy transfer ko = AE/c 
is of the same order as the momentum transfer k, the parameter is rather a/n. 
In the QED loops, such as for the vacuum polarization and the self-energy, the 
relativistic condition k ~ |k| is always a case. Sometimes it is also realized 
for the exchange loop and in particular in the case of some recoil corrections. 
However, most of the Coulomb exchanges correspond to another kinematics, 
which is close to the external field case with a low energy transfer ko <C |k| and 
a real parameter is even not Za but rather nZa at least for a few first terms 
of expansion. One of sources of the factor 7r is a fact that the imaginary part 
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of logarithm is typically n, not unity, and with multilogarithmic contributions 
some factors 7r in non-logarithmic terms must be expected. 

Because of a presence of at least three small parameters it is not possible to do 
any exact calculations and one must use expansions at least in some parame- 
ters. The Za expansion with all its enhancements has a very bad convergency 
and only in a few cases (some one-loop QED contributions [62] and recoil [63]) 
calculations exact in Za have been performed. In such the hardest the- 

oretical problem is to estimate systematically all uncalculated terms related 
to the higher-order corrections of the Za expansion. 

References to recent theoretical results on the Is Lamb shift theory are col- 
lected in Table 3. Calculation of these orders has been completed. Some other 
corrections are under study (e.g., a 3 (Za) 5 mc 2 [78]). 

We note that accuracy of the most recent calculations of the two-loop self- 
energy [67,68] is unclear and we ignore it in our analysis. In particular we 
doubt that the higher order terms may neglected in extrapolation to low Z 
[68] . On contrary we expect they can produce corrections bigger than claimed 
in [68] uncertainty of the calculation of the non-logarithmic term. We also do 
not think that a reference to a similarity with the one-loop correction is enough 
to claim a 15% uncertainty for the non- logarithmic term in [68]. The one-loop 
and two-loop self energies are quite different and in particular they show a very 
different n dependence when comparing contributions for various ns states. 
The n dependence is strongly related to splitting of the complete result to 
partial contributions of specific areas of integration. Besides we were recently 
informed that the single- logarithm two-loop calculation needs corrections [77]. 
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Table 3 

References to recent theoretical results on calculations of the Is Lamb shift. 



A pure QED calculation for the hydrogen atom is not enough and before 
performing any comparison of the QED theory and experiment we need to 
take into account effects of the internal structure of a proton. 
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5 Lamb shift and effects of the proton charge distribution 

Let us consider a problem related to the distribution of the proton charge in 
more detail. The leading correction due to the finite size of the proton has a 
simple expression 

, 2(Za) 4 mc 2 fmcR p \ 2 . 
A£ finitc sizc (nZ) = -^-^ {~T^J 6 ' 

where R p is the mean-squared proton charge radius. However, to reach any 
accurate numerical value one needs first to determine the value for the pro- 
ton charge radius with a proper precision. Unfortunately, data available at 
the moment, if understood literally with their claimed uncertainties, are not 
reliable (we follow here a consideration in [36]). The dominant contribution 
to the uncertainty of theory of the Lamb shift of the Is state is due to the 
contribution (19) [36]. It appears that currently the most accurate value of 
R p can be obtained by reading a comparison of the experimental value of the 
Lamb shift and theory in a reversed way: not testing the bound state QED, 
but determining the proton charge radius [36] 

i? p (Lamb shift) ~ 0.89(2) fm . (20) 

The prediction for the theoretical value of the Lamb shift strongly depends on 
what value for the proton charge radius is accepted. A crucial level of accuracy 
of the radius is 1%, which is related to the case when all three uncertainties 
(experimental, QED and due to the proton radius) are approximately the 
same. Two values were under a long discussion inside the atomic and QED 
community: 

i? p (Stanford) = 0.805(11) fm , [79] , (21) 

and 

i? p (Mainz) = 0.862(12) fm , [80] 

More data for the proton radius are collected in Fig. 13. 

To discuss the scatter of the data, let us look at the most important data on 
the electron-proton elastic scattering presented in Fig. 14. Even from the first 
glance one can realize that the Mainz data are 

• presented by a larger number of points; 

• more accurate; 
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Fig. 13. Proton charge radius determined from the scattering experiments. The 
presented results are phenomenologically extracted from the scattering at Orsay 
[81], Stanford [79], Saskatoon [82] and Mainz [80] or found with more sophisticated 
analysis from Mainz data by Wong [83] and by Mainz theoretical group from a 
multi-parameter dispersion fit of all available data [84]. 

• obtained for a more broad range of the momentum transfer including the 
lowest momentum transfer ever achieved. 



The latter is important since the easiest way to extract the proton charge 
radius from the data is to extrapolate to zero momentum transfer 
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Fig. 14. Data for the electric form factor of the proton from the electron-proton 
elastic scattering experiments performed at Orsay [81], Stanford [79], Saskatoon 
[82] and Mainz [80]. 
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where q 2 = — q 2 is related to the three-dimensional momentum transfer. There 
is no chance for other data in Fig. 14 to compete with the Mainz results [80] 3 . 
The Mainz experiment is the most appropriate to precisely determine the 
proton radius since it contains more points at lower momentum transfer and 
with a higher precision than the rest of the data. Due to this any compilation 
containing the Mainz data has to lead to a result close to that from the Mainz 
data only, because the Mainz scattering points must be statistically responsible 
for the final result. In particular, we believe the dispersion analysis performed 
by Mainz theorists [84] led to such a result 4 

i2p(disp) = 0.847(9) fm , (24) 



which, however, differs from the empirical value of R p = 0.862(12) fm from 
Eq. (22) [80]. 
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Fig. 15. Fitting the electric form factor of the proton from the Mainz experimen- 
tal data [80]. Since details of fitting by Wong [83] and a particular result on the 
coefficient 02 are not available, we present here similar fits from [36]. 

One of the problems in evaluating the data is their normalization and actually 
a different treatment of this problem is in part responsible for the discrepancy 
between results from Eqs. (22) and (24). One can write a low- momentum ex- 
pansion of the form factor from Eq. (23) and expect that G(0) = 1. However, 
a real situation is more complicated and somewhat confusing. The confusion 
arises from use of the words 'form factor'. The definition of the form factor is 

3 I have learned that from a similar picture presented by Malcolm Boshier for ICAP 
1996 (see Fig. 5 in [85]). 

4 Unfortunately, the authors of [84] are neither able to explain which part of data 
statistically dominates (if any) nor to prove properly the estimation of the uncer- 
tainty assigned by them to the nucleon radii since the radii were not main objectives 
of their fitting procedure and appeared in the very end of the evaluation as a com- 
plicated function of numerous correlated fitting parameters. 
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clearly understood theoretically, however, the true form factor cannot be mea- 
sured accurately. In the case of scattering one has to deal with the Rosenbluth 
formula which describes the cross section. To make an absolute measurement 
of the cross section, one needs to measure parameters of the incoming electron 
beam with a high accuracy, as well as determine the efficiency of the detectors 
etc. For a relative measurement one has only to ensure that the details of the 
beam and detections are the same during the entire experiment. The normal- 
ization is very complicated and sometimes involves much more efforts than 
the relative measurements. However, the normalization measurement used to 
be not enough accurate. In particular, in the Mainz case the uncertainty is 
about 0.5%. That means that a value of G(q 2 ) tabulated from the scattering 
data, as being the form factor, actually differs from the true form factor with 
some normalization factor, which is consistent with but not equal to unity. 

There is an alternative wording to explain the situation and it may be helpful 
for some readers. When the experimentally determined form factor is pre- 
sented properly it contains a statistical error (mainly related to relative mea- 
surements) and a systematic error (mainly due to the normalization). The 
measured form factor in such a case directly related to the 'theoretical form 
factor' and its value at zero momentum transfer is equal to unity. However, 
as we mentioned the systematic error is often higher than the statistical er- 
ror (at least eventually, when the statistical error for the fitting parameters 
is suppressed by the number of data points). That means that if we at the 
first approximation would ignore the systematic error we still should be able 
to obtain a result almost without loosing in accurate. If we drop the system- 
atic error, the central values of the form factor remain the same, however, it 
physical meaning is changed drastically. The measured value G meas (q) can be 
presented as 

G mcas (q)=G' mcas (q)-(l±5), 

where 5 is the systematic error of the normalization. While the measured 
form factor is defined in such a way that G meas (0) = 1; the value with only a 
statistical uncertainly satisfies a condition G' meas (0) = 1/(1 ± 5), i.e., it must 
be consistent with unity at zero momentum transfer, but not just equal to it. 
That is to what we refer above as to a conventional normalization. 

The confusion comes in part from the common attempt to normalize data as 
close to the form factor as possible to make their physical meaning more clear. 
However, from a pure practical point of view that has no sense. The lack of ac- 
curate data on the normalization makes G(0) consistent with unity within the 
experimental uncertainty of the normalization, which is nearly always larger 
than the uncertainty of relative measurements. The normalization uncertainty 
is a systematic error, while the data are often tabulated with only statistical 
component of the uncertainty. 
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A few different fits were performed by Wong [83] (see Fig. 15). Two of them 
were related to the normalization used in [80] and [84] 5 and the achieved 
results were close to published there. A fit with a free value of a led to a 
larger uncertainty R p = 0.88(2) fm (see the Wong-Mainz value in Fig. 13). 
To our mind, even this result must be treated with a certain caution. It is 
necessary to take into account some higher-order corrections and that is not 
possible because of the absence of any complete description of the experiment 
[80]. Our conservative estimation of the theoretical uncertainty in Fig. 9 is 
related to [36,60] 

^(conservative) = 0.88(3) fm . (25) 

The present status of the Lamb shift of the Is (2s) state is that 

• the QED computation uncertainty is about 2 ppm, 

• the measurement uncertainty of the grand average value is 3 ppm, 

• while the uncertainty due to the proton size is about 10 ppm. 

Some other analysis of the proton radius can be found in [86]. The result 
R p = 0.880(15) fm is consistent with (25). It is likely that a proper normaliza- 
tion still is the most important source of correcting the central value of R p and 
the increase of its uncertainty. A very different treatment by Sick [87] which 
takes into account Coulomb effects and various systematic sources led to a 
result R p = 0.895(18) fm also consistent with our value in (25). Our uncer- 
tainty is somewhat bigger than that of [83,86,87] because we consider possible 
systematic errors due to higher order QED corrections [36] to scattering which 
are beyond evaluation in [80]. 

Improvement of the determination of the proton charge radius is expected from 
a study of the Lamb shift in muonic hydrogen. As one can see from discussion 
in Sect. 2, the line width scales linearly with the mass of the orbiting particle, 
while the finite-nuclear-size contribution is proportional to m 3 and that offers 
a good chance for an accurate determination of the proton charge radius. The 
level scheme of the experiment which is in progress at Paul Scherrer Institut 
[88] is presented in Fig. 16. The metastable 2s state is to be excited by the 
laser light to the 2p state which then decays to the ground state producing an 
X-ray Lyman-o; photon. The intensity of the X-ray fluorescence is studied as 
a function of the laser frequency. 

The problem of the nuclear size is not only a problem of the hydrogen Lamb 
shift: a similar situation arises with the helium-4 ion Lamb shift, where uncer- 
tainties resulting from the QED computation and the nuclear size are about 

5 The paper [83] was published prior to [84], however, it was indeed natural to 
anticipate the use of the 'theoretical' normalization G(0) = 1. 
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Fig. 16. Level scheme of the PSI experiment on the Lamb shift in a muonic hydrogen 
[88] (not to scale). The hyperfine structure is not shown. 

the same. Theory is in agreement with experiment [89] and the leading theo- 
retical uncertainty is related to the nuclear size. The theoretical prediction for 
the Lamb shift in deuterium has a bigger uncertainty because of the nuclear 
polarizability effects [90]. The latter are significant since a deuteron is known 
as a very loosely nuclear bound system. For more detail on the nuclear physics 
of light nuclei and their effects on atomic energy levels, see [91]. 



6 Hyperfme splitting in light hydrogen-like atoms and the nuclear 
structure 

A similar problem of significant interference of the nuclear structure and QED 
effects exists for the hyperfme structure (HFS) in the hydrogen atom and ac- 
tually the case of the hyperfine structure is in a sense even much less favorable 
for the QED tests. The hyperfine interval in the ground state was for a while 
the most accurately measured quantity being known with an uncertainty at 
the level of a part in 10 12 . The results of the most accurate experiments are 
summarized in Fig. 17. It so happened that most results were published in 
metrological journals and were missing in various 'physical' compilations. This 
figure includes only results published in refereed journals, but not those from 
conference proceedings. We also note that the most frequently quoted result 
[102] 

/hfs(Is) = 1420 405.751767(1) kHz 

was later corrected [94] and the uncertainty was increased by a factor of 3. 

We summarize in Tables 4 all accurate data available for the hyperfine intervals 
of the Is and 2s states in conventional light hydrogen-like atoms (hydrogen, 
deuterium and tritium and in the helium-3 ion). Pure leptonic atoms are 
considered separately in Sects. 8 and 9. 
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Fig. 17. Measurement of the Is hyperfine splitting in the hydrogen atom. The ref- 
erences are: a - [92] (experiment #1 with the wall-shift determined in [93]), b - [92] 
(experiment #2), c - [94], d - [95], e - [96], / - [97], g - [98], h - [99], i - [100]. The 
conservative average, 1420 405.751 768(1) kHz, is presented according to [101]. It is 
an average over the results denoted by filled circles, while the results presented by 
open circles are excluded from evaluation (see [101] for detail). 

A value of the hyperfine interval in the ns state is determined in the leading 
order by the so-called Fermi energy Ep as following 



E ns (leading) = E F 



n 3 



8na <Me-/% ucl > |vMr = 0)|2 _ (26) 



3 mlc ii 2 B 

Here fi e and /iNuci stand for the magnetic moment of a free electron and the 
nucleus, respectively, and 



is a squared value of the non-relativistic wave function at the origin. Note, 
that we have included a contribution of the anomalous magnetic moment of 
the electron into the leading term E ns (leading) , but not into the definition of 
Ep. 

The most important corrections to (26) are due to relativistic effects, bound- 
state quantum electrodynamics and due to the nuclear structure. In the case of 
a point-like infinitely heavy nucleus, i.e., when only the relativistic and bound 
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Atom, 


-^HFs(exp) 


Ref. 


Ehfs(QED) 


AE 


state 


[kHz] 




[kHz] 


[ppm] 


Hydrogen, Is 


1420 405.751768(1) 


[92,93,94,95,98,99,101] 


1420 452 


- 33 


Deuterium, Is 


327 384.352 522(2) 


[103] 


327339 


138 


Tritium, Is 


1516 701.470 773(8) 


[104] 


1 516 760 


- 38 


3 He + ion, Is 


- 8 665 649.867(10) 


[105] 


- 8 667494 


- 213 


Hydrogen, 2s 


177556.860(16) 


[106] 


177 562.7 


- 32 


Hydrogen, 2s 


177 556.785(29) 


[107] 




- 33 


Hydrogen, 2s 


177 556.860(50) 


[108] 




- 32 


Deuterium, 2s 


40 924.454(7) 


[109] 


40 918.82 


137 


Deuterium, 2s 


40924.439(20) 


[110] 




137 


3 He + ion, 2s 


- 1083 354.980 7(88) 


[111] 


- 1083 585.7 


- 213 


3 He+ ion, 2s 


- 1083 354.99(20) 


[112] 




- 213 



Table 4 



Hyperfine structure in light hydrogen-like atoms. The numerical results are pre- 
sented for the frequency E/h. The difference (AE) between the actual value 
(-£'HFs( ex P)) an d a result of the pure QED calculation (E'hfs(QED)) is caused by 
the nuclear structure. The negative sign for the 3 He + ion reflects the fact that the 
nuclear magnetic moment is negative, i.e., in contrast to other nuclei in the Table, 
its direction is antiparallel to the nuclear spin. 

state QED contributions are taken into account, the correction to (26) is of 
the form 



A£ ls (QED) = Ep x \^{Zaf + a(Za) (hi 2 
a(Za) 2 



+ 

7T 

+ 41n2 

8 



In 



3 (Zaf \ (Zaf 



281 
240 
34 



+ 17.122 339. . . 



In 2 + 

15 225 



+ 0.7718(4)^1} 



(28) 



This term is in fact smaller than the nuclear-structure correction (see Table 5). 
The contributions induced by the nuclear structure are determined there by 
subtracting the results of the pure QED calculations 



Ehfs(QED) = £ ns (leading) + A£ ns (QED) 



41 



from the actual values of the hyperfine intervals (see Table 4) 



Ehfs(NucI) = £ H Fs(exp) - £ H fs(QED) . 



The result for the Is state is presented in Eq. (28) and the 2s result is of a 
similar form and discussed in detail in the next section. We mainly follow here 
[113], but the QED helium result in Table 4 and related data in Tables 5 and 
6 are corrected as explained in [114] and [115,116] (cf. [113]). 



Atom 


AE'(QED) 


{Za) 2 E F 


a(Za)Ep 


A£(Nucl) 




[ppm] 


[ppm] 


[ppm] 


[ppm] 


Hydrogen 


-23 


80 


-96 


-33 


Deuterium 


-23 


80 


-96 


138 


Tritium 


-23 


80 


-96 


-38 


3 He + ion 


108 


319 


-192 


-213 



Table 5 

A comparison of the bound QED and nuclear-structure corrections to the Is hyper- 
fine interval. The QED term AE'(QED) defined in (28) contains only the bound- 
electron effects while the anomalous magnetic of free electron a e = a/2ir + ... ~ 
1.159 652 ... x 10~ 3 is excluded. Additionally, to the whole QED correction we also 
present separately two major contributions into it: the leading relativistic term in 
order (Za) 2 Ep and the leading bound state QED correction in order a(Za)Ep- The 
nuclear contribution Ai£(Nucl) has been found via a comparison of the experimental 
results with the pure QED values (see Table 4). 



There are three kinds of corrections which involve the nuclear-structure effects. 



The dominant contribution for the hydrogen HFS interval is due to a dis- 
tribution of the charge and magnetic moment inside the nucleus (so-called 

Zemach correction) 



Ai^(Zemach) = — ^- E F x mc [ 
it 2 J 



d 3 q 



G E (-q 2 )G M (-q 2 ) 

1 + K 



- 1 



, (29) 



where the electric/magnetic form factor Ge/m(— Q 2 ) is a Fourier transform 
of the space distribution of the proton electric charge/magnetic moment 
Pe/m(j) and k is the proton anomalous magnetic moment 6 . In the coordi- 



6 It is necessary to mention that a straightforward consideration of the form fac- 
tor as a Fourier transform of the space distribution of the related value (such as 
the electric charge) is a transparent but rough approximation, which fails for the 
relativistic momentum transfer (see, e.g., [117]). 
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nate space the equation takes the form 

AE(Zemach) = _ 2( ^ Zo ^ mc E F x J dhdh' p E (r) PM (r') |r - r'| . (30) 

For light atoms the uncertainty in (29) comes from the distribution of the 
magnetic moment which involves the magnetic form factor Gm(— q 2 ) at 
low momentum transfer, which is difficult to study experimentally in such 
kinematic region (in contrast to the electric form factor Ge)- 

• A different effect dominates for a deuteron, which is a loosely bound system 
and the electron can see it not only as the whole but can also recognize its 
composites, namely a proton and a neutron. This is a nuclear polarizability 
contribution. 

• A theory of a point-like particle with a non- vanishing value of the anomalous 
magnetic moment is inconsistent 7 . In particular, calculating recoil correc- 
tions to the hyperfme structure in atoms with such a nucleus, certain diver- 
gencies have to appear at high momentum. Those are not ultraviolet QED 
divergencies. They are related to an inconsistency of the consideration of the 
operators of the non-relativistic expansion as fundamental operators. The 
operator a^k^/Mc is such an operator and it gives reasonable results only 
when k/Mc < 1. A complete theory with the finite nuclear size removes 
all divergencies of this kind and thus the nuclear size serves as an effective 
cut-off for them entering the equations originally designed for a point-like 
particle. Since the divergence in the related recoil contribution to the HFS 
interval is logarithmic, the recoil effects only slightly depend on this cut-off 
and through it on detail of the nuclear structure. 

The magnitude of the nuclear corrections to the hyperfme interval in light 
hydrogen-like atoms lies at the level from 30 to 200 ppm (depending on the 
atom) and theoretical understanding of such effects is unfortunately by far 
not sufficient [36,118,113]. The correction is compatible with the bound-state 
QED term (note that the anomalous magnetic moment is related to the free 
QED consideration and it is included into the leading term, but not into the 
QED correction). The results for the Is state are summarized in Table 5 (see 

7 We need to mention a certain inconsistency in the terminology related to the 
leptons. They are referred as "point-like" since all effects of the distribution of the 
electric charge and the magnetic moment can be derived from QED. However, they 
"have a non-vanishing anomalous magnetic moment" which comes from the same 
QED effects. Indeed, both effects are to be considered on the same ground. The in- 
consistency historically comes from the fact that the anomalous magnetic moments 
was first discovered experimentally and later explained, while the "electron charge 
radius" is hardly to be seen in a straightforward experiment but can be easily con- 
sidered theoretically. The Lamb shift is a result of the electron charge distribution 
related to the electron structure, however, it is understood only through the QED 
interpretation. So the internal electron structure effects look rather as a theoretical 
construction, while the anomalous magnetic moment is a direct experimental fact. 
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Ref. [113] for detail). The results for the 2s state are slightly different for the 
relative value of the QED contribution. A certain progress [119] is on the way, 
however, it should more contribute to nuclear physics since a gap between 
accuracy of a pure QED theory and that of the nuclear effects is of a few 
orders of magnitude and cannot be avoided. 

From Table 4 one can learn that in the relative units the effects of the nuclear 
structure are nearly the same for the Is and 2s intervals (33 ppm for hydrogen, 
138 ppm for deuterium and 213 ppm for helium-3 ion). A reason for that is 
the factorized form of the nuclear contributions in the leading approximation 
(cf. (19)) 

AE(Nucl) = A(Nucl) x |tf nZ (r = 0)| 2 . (31) 



In other words, the correction is a product of the nuclear-structure parameter 
A(Nucl) and the squared wave function at the origin (27), which is a result 
of a pure atomic problem (a nonrelativistic electron bound by the Coulomb 
field). The nuclear parameter A(Nucl) depends on the nucleus (a proton, a 
deuteron etc.) and on the effect under study (the hyperfine structure, the 
Lamb shift), but does not depend on the atomic state. We remind, that, 
similarly, the leading term for the HFS interval for an ns state in Eq. (26) 
is also proportional to squared wave function at the origin ^^(r = 0)| 2 and 
thus a fractional value of the nuclear contribution is kept n independent in 
the leading approximation. 

Two parameters in the hydrogen-like wave function (27) can be varied for the 
same nuclei: 

• the reduced mass of a bound particle for conventional (electronic) atoms 
(rriR ~ m e ) and muonic atoms (m,R ~ m M ); 

• the principle quantum number n — 1, 2 for the Is and 2s states. 

We note that the nuclear charge Z is indeed fixed if we like to deal with the 
same parameter A(Nucl) (i.e., with the same nucleus). 

The former option was already mentioned when we considered the determi- 
nation of the proton charge radius via the measurement of the Lamb shift 
in muonic hydrogen [88]. In the next section we consider the latter option, a 
comparison of the Is and 2s hyperfine intervals in hydrogen, deuterium and 
the 3 He+ ion. The Is state is the ground state and the 2s state is metastable. 
Other levels are not accessible for the precision experiments being broad be- 
cause of their radiative decay. 
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7 Hyperfine structure, free of the nuclear effects: comparison of 
the HFS intervals for the Is and 2s states 



As we have seen in the previous section, both the Is and 2s hyperfine intervals 
in various light hydrogen-like atoms are measured more accurately than theory 
can predict. Theory, although limited by the uncertainty in the calculation of 
the nuclear effects, can still proceed to a higher accuracy dealing with a specific 
difference 

D 21 = 8 ■ E HFS (2s) - E HFS (ls) , (32) 



where any contribution (and, in particular, related to the leading nuclear ef- 
fects) which has the form of Eq. (31) should vanish. Removing the main ob- 
stacle on the way of theory, we may study effects related to higher-order QED 
contributions. 

The difference (32) was studied theoretically in several papers long time ago 
[120,121,122,123]. A result of the former study reads 



D 21 (QED3) = {Zaf E F x ( - + - [ ( — In 2 - 7^ ln(Za) - 5.551 55 . . . 

18 % LV 3 / 



+ 



m 
M 



where 



li M 1 

V = 1TJ 

Hb m Z 1 



ln2 



32 



16 
3~ 



145 7 



128 8 



In 2 ) 77 



(33) 



(34) 



We remind here about a difference in the notation in particle and nuclear 
physics. If the g factor of a nucleus would be defined in the particle notation, 
one has to arrive at r\. 

The HFS difference D 2 \ is similar to A(n) in Eq. 11, however, there are a num- 
ber of very different features. The HFS difference is designed to be compared 
with experimental data and it is of interest only for the Is and 2s states. The 
Lamb difference A(n) is to be used for the interpretation of the data and those 
involves various n, still n = 2 is the most important for application case. The 
leading term E ns {leading) , the leading free QED (due to the anomalous mag- 
netic moment of the electron) and the leading bound state QED corrections 
cancel for the HFS difference. On the contrary, the leading contribution to the 
Lamb shift does not cancel for A(n). Still, there is a similarity in calculations 
of various bound state QED contributions. 
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L> 2 i(QED4) 


0.018(5) 


0.004 4(10) 


-1.13(14) 


L> 21 (Nucl) 


-0.002 


0.002 6(2) 


0.307(35) 


D 2 i(theo) 


48.953(5) 


11.312 5(10) 


-1190.08(15) 



Table 6 

Theory of the specific difference D21 = 8.Ehfs(2s) — -Ehfs(Is) in light hydrogen-like 
atoms (see [115,116,113,114] for detail). The numerical results are presented for the 
related frequency D21/I1. 

Recent studies of D21 [124] showed that certain higher-order QED and nu- 
clear corrections have to be taken into account for a proper comparison of 
theory and experiment. The theory has been essentially improved [113,125] 
and it is summarized in Table 6. The QED corrections up to the third order 
(-D2i(QED3)) and the fourth-order contribution of the order (-Za) 4 (in units of 
the Is hyperfine interval) have been known for a while [121,122,123,120]. The 
new issue here is a complete result for the fourth-order QED contributions 
(_D 2 i(QED4)) and nuclear corrections (_D 2 i(nucl)). 

The QED contributions in the fourth order are summarized in Table 7. The 
new QED contributions are of the order a(Za) 3 , a 2 (Za) 4 , a(Za) 2 m/M and 
(Za) 3 m/M. The difference with [113] is a more realistic estimation of the 
a(Za) 3 uncertainty [115,116]. In the former paper we accepted result from 
[125]. 

Both results for a(Za) 3 at Z — 1,2, in [115,116] and [125], are based on 
extrapolations of numerical data obtained in [125] for higher values of Z. The 
procedure leads to a calculation of central values and its uncertainty. While 
the calculation of the central value is approximately the same, the uncertainty 
is quite different. Systematic uncertainty of the fitting pro cedure is due to a 
possible perturbation of the fit by certain higher-order terms. As we shown in 
[115,116], one should expect certain relatively large numerical coefficients for 
these terms and that increases the uncertainty. 

The higher-order nuclear effects survive the cancellation of the leading term 
but they are substantially smaller than the leading term and can be success- 
fully calculated phenomenologically [113] 



fl2 1 =(ln2 + -).(Z a ) 2 .£ h N ,r , (l S ) 

+ (2 - £ ln 2 ) • (Z„)» i^)E F - 1 • (ZaY (^fif E F .(35) 
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Contribution 


H 


D 


3 He+ 


Ref. 




[kHz] 


[kHz] 


[kHz] 




{Za) 4 E F 


0.005 6 


0.0013 


-0.543 


[120] 


a 2 {Za) 2 E F 


0.003 3(16) 


0.000 8(4) 


-0.069(35) 


[124,113] 


a(Za) 2 ^E F 


-0.0031(15) 


-0.000 4(2) 


0.022(11) 


[124,113] 


a(Za) 3 E F (SE) 


0.008(4) 


0.0019(9) 


-0.39(14) 


[115,116,125] 


a{ZafE F (VP) 


0.003 


0.000 7 


-0.145 


[124,113] 


(Zaf§E F 


0.000 5(5) 


0.000 1 


-0.007(10) 


[124,113] 


£» 2 i(QED4) 


0.018(5) 


0.004 4(10) 


-1.13(14) 





Table 7 

The fourth order QED contributions to the D21 in hydrogen, deuterium and helium- 
3 ion (see [115,116,113] for detail). 

The results for the light atoms are summarized in Table 6. A value for the 
nuclear correction in the helium-3 ion is slightly different here from the result 
in [113] due to different nuclear parameters we used (see Table 8). 

Atom Z I M/m m/M-B E f R f Rm 

[10~ 3 ] [kHz] [fm] [fm] 

Hydrogen 1 1/2 1836.153 1.521032 1418 840 0.88(3) 0.86(4) 

Deuterium 1 1 3 670.483 0.466 9755 326 968 2.13(1) 2.07(2) 

Tritium 1 1/2 5 496.922 1.622 3936 1515 038 1.76(9) 1.84(18) 

3 He+ 2 1/2 5 495.885 -1.158 741 -8 656 527 1.96(3) 1.97(15) 

Tabled 

Parameters for the calculations of the HFS interval in hydrogen, deuterium and 
helium-3 ion. The proton charge radius is presented according to Ref. [36] while the 
magnetic radius is from [126]. The radii for other nuclei are taken from [127], while 
the masses and magnetic moments are from [23,128,129,115]. One should use the 
values for the radii with caution. The results in [127,87,126] are obtained from the 
scattering data as they were published. However, in a number of papers published 
decades ago certain radiative corrections were not taken into account (cf. our analy- 
sis in [36] ) and even a perfect evaluation with uncorrected data underestimates the 
uncertainty. 

For all atoms in Table 4 the hyperfine splitting in the ground state was mea- 
sured more accurately than for the 2s state. All experimental results but two 
were obtained by a direct measurement of the microwave transitions for the Is 
and 2s hyperfine intervals. However, two of the most recent results obtained 
for hydrogen [106,115] and deuterium [109,115] have been achieved by means 
of laser spectroscopy and the measured transitions lie in the ultraviolet range 
[106]. The hydrogen level scheme is depicted in Fig. 18. The measured transi- 
tions were the singlet-singlet (F = 0) and triplet-triplet (F = 1) two-photon 
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2S1/2- 



2shfs(rf) 



two-photon 
uv transitions 



1S 



1/2 



1s hfs (rf) 
t 



F = 1 (triplet) 



F = (singlet) 



Fig. 18. Level scheme for the optical measurement of the 2s hyperfine structure in the 
hydrogen atom (not to scale) [106,115]. The label rf stands here for radiofrequency 
intervals, while uv is for the ultraviolet transitions. 

Is — 2s ultraviolet transitions. The eventual uncertainty for the hyperfine 
structure is related to 5 parts in 10 15 of the measured Is — 2s interval. A sim- 
ilar experiment for deuterium achieved an accuracy of the hyperfine interval 
related to even a smaller portion of the 'big' Is — 2s transition, namely 3 parts 
in 10 15 . 
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Fig. 19. Present status of experiment and theory of D21 in the hydrogen atom. The 
results are labelled with the date of the measurement of the 2s hyperfine structure. 
See Table 4 for references. 

The comparison of theory to experiment for hydrogen, deuterium and helium- 
3 ion is summarized in Figs. 19, 20 and 21. We acknowledge a substantial 
recent progress in experimental studies of neutral hydrogen and deuterium 
atoms, however, the helium ion data are still of the most interest since theory 
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Fig. 20. Present status of experiment and theory of D21 in the deuterium atom. The 
results are labelled with the date of the measurement of the 2s hyperfine structure. 
See Table 4 for references. 

and experiment are compatible. 
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Fig. 21. Present status of experiment and theory of D21 in the helium ion 3 He + . The 
results are labelled with the date of the measurement of the 2s hyperfine structure. 
See Table 4 for references. 

8 Hyperfine splitting in the ground state of muonium 

Another possibility to eliminate the nuclear structure effects is based on stud- 
ies of nucleon-free atoms. Such an atomic system can be formed of two leptons. 
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Two atoms of the sort have been produced and studied for a while with high 
accuracy, namely, muonium and positronium. We consider the former below 
in this section, while the next section is devoted to the latter. 

Muonium is a bound system of a positive muon and an electron. It can be 
produced with the help of accelerators. The muon lifetime is 2.2 • 10 -6 sec. 
However, it is sufficient to perform accurate spectroscopic measurements. The 
most accurately measured transition is the two-photon Is — 2s transition with 
an uncertainty of 4 ppb [130]. The other accurately measured interval is the 
Is hyperfine splitting (12 ppb) [131]. A detailed review of muonium physics 
can be found in [132]. 



Term 


Fractional 


AE 




contribution 


[kHz] 


Ep 


1.000 000 000 


4459 031.88(50)(3) 


a e 


0.001 159 652 


5170.925(1) 


QED2 


- 0.000 195 815 


- 873.145 


QED3 


- 0.000 005 923 


- 26.411 


QED4 


- 0.000 000 123(49) 


- 0.548(218) 


Hadronic 


0.000 000 054(1) 


0.241(4) 


Weak 


- 0.000 000 015 


- 0.067 


Total 


1.000 957 830(49) 


4 463 302.88(51)(3)(22) 



Table 9 

Theory of the Is hyperfine splitting in muonium. The calculations [133] have been 
performed for a' 1 = 137.035 998 76(52) [134] and av/^p = 3.183 34517(36) which 
was obtained from the analysis of the data on the Breit-Rabi levels in muonium 
[131,135] (see Sect. 10.4) and muon precession [136]. The numerical results are 
presented for the frequency E/h. 

Despite that the fractional accuracy of the measurement of the hyperfine in- 
terval is three times smaller than that of the Is — 2s transition, it is of a 
bigger interest. The uncertainty of the hyperfine measurement in muonium is 
indeed much smaller in absolute units (53 Hz against 9.8 MHz for the Is — 2s 
interval) and thus more sensitive to higher-order corrections. The Lamb shift 
of the Is state in muonium can be extracted from the optical measurement of 
the Is — 2s transition with fractional uncertainty of approximately 1%, which 
is indeed not enough for any interesting QED test (it is larger, e.g., than the 
proton-finite-size contribution in the case of hydrogen). A real importance of 
the Is — 2s experiment is due to the determination of the muon mass which 
can be used in theoretical calculations for other QED tests such as the muo- 
nium hyperfine structure. We discuss the muon mass and magnetic moment 
in Sect. 10.4. 
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Let us consider now the hyperfine interval of the ground state in muonium. 
The theoretical status is presented in Table 9. The budget of theoretical un- 
certainties in the case of muonium involves quite different sources and they 
are in part of experimental origin. The biggest source is in the calculation of 
the Fermi energy, accuracy of which is limited by the knowledge of the muon 
magnetic moment or the muon mass in proper units. To determine the muon 
magnetic moment (<7 M /2) (eh/2m fl ) and its mass is essentially the same because 
the g factor of a free muon is known well enough [137]. A smaller uncertainty 
in determination of the Fermi energy (see the number in the second brackets) 
is due to the fine structure constant. The theory up to the third order (i.e., 
a e , QED2 and QED3 terms) is well known (see references in Table 10). The 
nonrecoil part up to the third order is determined for the ground state by (26) 
and (28), while the recoil effects lead to an additional contribution 
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Correction 


References 


Correction 


Reference(s) 


a(Za) 2 E F 


[138,139,140] 


(Za) 2 (m/M)E F 


[142] 


a 2 (Za)E F 


[141,66] 


a{Za)(m/M)E F 


[143,144,145] 



Table 10 

Recent references on the third order QED corrections for the ground state muonium 
hyperfine splitting. See [16] for detail. 

The fourth-order corrections (AE(QED4)) have not yet been included either 
into a non-recoil term E\ S (QED) (see (26) and (28)) or into a recoil term 
AEi s (recoil). They are the corrections that determine the uncertainty related 
to QED theory. These corrections are related to essentially the same diagrams 
as the D2\(QEDA) contribution in the previous section. The results are sum- 
marized in Table 11. Unknown corrections are estimated according to [133]. 
Meantime some of them and in particular the corrections of relative order 
a 3 (Za) [78] and a 2 (Za)(m/M) [146] are under study and a few gauge in- 
variant sets of diagrams have been already calculated. The muonium QED 
uncertainty is mainly due to the calculation of the recoil corrections of order 
a(Za) 2 m/M and (Zafm/M [148]. 

The muonium theory is not completely free of hadronic contributions, which 
appear because of intermediate hadronic states. That produces one more 
source of theoretical uncertainties. The effects of strong interactions in muo- 
nium are discussed in detail in Refs. [133,161,162]. They are relatively small, 
but their understanding is very important because of the intensive muon 



51 



Contribution Numerical result 



Reference 



a{Zaf -0.52(3) kHz 

a(Za) 2 § 0.39(17) kHz 

a(Za) (f f -0.04 kHz 

a 2 (Za) 2 -0.04(2) kHz 

a 2 (Za)^ -0.04(3) kHz 

a 3 (Za) ±0.01 kHz 

r lepton 0.002 kHz 




-0.29(13) kHz 
-0.02(1) kHz 



0.03 kHz 



0.02 kHz 



[154] 

[155,125,156], [148], 



[61,147,141,148,149 
[151,152,153] 



[120] 



150] 



[157] 



[61,149,150] 



[158] 
[61] 
[159,160] 

[133] 
[161,162] 



A£(QED4) -0.55(22) kHz 



Table 11 



The fourth order corrections to the muonium hyperfine structure. The uncertainty 
is an rms sum of the partial uncertainties. The essential part of the contributions 
(see [133] for detail) have been found within the leading logarithmic approximation 
and their uncertainty is estimated for each contribution as a half- value of the leading 
term [148,133]. The r lepton contribution is also included into the table and has an 
order a(Za)mM/m 2 . 

sources expected in future [163] which can allow a dramatic increase of the 
accuracy of muonium experiments. 

A comparison of theory and experiment for the muonium HFS interval is not 
an isolated problem, since it involves determination of the muon mass (mag- 
netic moment) in proper units. We discuss problems of the determination of 
the muon mass and magnetic moment in Sect. 10.4 and a comparison of theory 
versus experiment for muonium hyperfine interval is presented in the Sect. 13 
along with the other test of bound state QED for the hyperfine structure (see 
Table 19). 



9 Spectrum and annihilation of positronium and recoil effects 

Another pure leptonic atom is positronium. It can be produced at accelerators 
or using various radioactive positron sources. The lifetime of positronium de- 
pends on its state, on its orbital and spin (S) quantum numbers. The lifetime 
for the IS state of parapositronium (this state with S = annihilates mainly 
into two photons) is 1.25 • 10~ 10 sec, while orthopositronium (the S — 1 state) 
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in the IS* state has a lifetime of 1.4 • 1CT 7 sec because of the three-photon 
decays. A list of accurately measured positronium quantities contains the IS* 
hyperfine splitting, the IS* — 2S interval, the 2S — 2P fine structure intervals, 
the lifetime of the IS state of para- and orthopositronium and several branch- 
ing ratios of their decays. We note that since the nuclear spin effects are not 
suppressed, positronium has a structure of energy levels (in respect to their 
spin and angular momentum) rather similar to a two-electron system (such 
as the neutral helium atom) and thence we use capital letter for its orbital 
momentum. 

Theoretical predictions for most of the values of interest and references to 
their comparison to experimental data are summarized in Table 12. 



Quantity 


Leading term 


Prediction 


Figure 


Au(lS-2S) 


| a 2 mc 2 


1233 607 222.2(6) MHz 


24 


Av HFS (lS) 


7 4 2 

j2 a mc 


203 391.7(6) MHz 


23 


Au(2 3 S! - 2 3 P ) 


5 4 2 

a mc 


18 498.25(9) MHz 


28 


Av(2 3 S! - 2 3 P 1 ) 


7 4 2 
-j^g2 CL T7~hC 


13 012.41(9) MHz 


28 


Au(2 3 Si - 2 3 P 2 ) 


y2Q Ci TTXC 


8 625.70(9) MHz 


28 


Av(2 3 S! - 2 1 P 1 ) 


14 2 

24 a mc 


11185.37(9) MHz 


28 


r(p-Ps) 


\a b mc 2 


7989.62(4) /is" 1 


26 


r(o-Ps) 


2(tt 2 -9) 6 2 

v 9n ' a mc 


7.039 96(2) ^s- 1 


25 


Br 47/27 (pPs) 


0.0278 a 2 


1.439(2) • 10~ 6 


27 



Table 12 

Theoretical predictions for positronium. Comparison to experimental data is pre- 
sented in figures quoted in the table. The details and references for experiment 
can be found in [166,167] while for theory in [168,169] and below in Table 16. The 
leading term for each value is presented in energy units. 

Fig. 22. Two-photon basic diagrams for the calculation of recoil contributions. 

In the case of the atoms with a 'heavy' nucleus there is a significant difference 
between contributions to the Lamb shift and the hyperfine splitting. The hy- 
perfine effects used to affect the Lamb shift and related intervals (such as the 
gross and fine structure) only slightly because the hyperfine splitting is of the 
order (Z a) A m 2 c 2 / M , while the Lamb shift is of the order a(Za) 4 mc 2 ln(l/Za). 
Since the effects due to the hyperfine shift are small, the higher order correc- 
tions to the HFS interval are not important for a study of the Lamb shift 
and related transitions. Entering into detail, we note that a crucial problem 
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for theory of the hyperfine structure is a calculation of the recoil corrections, 
while for the Lamb shift they are not so important. The difference originates 
from the fact that if the momentum of integration is significantly below Mc, 
the nuclear recoil corrections are rather of a kinematic origin similar (but much 
more complicated) to the reduced mass effects in the relativistic case (for the 
electron), since the nucleus can still be treated as a non-relativistic particle. If 
the momentum of integration is of the order of Mc, the diagram for the basic 
block (see Fig. 22) has order (Za) 5 (m/M) 2 mc 2 and higher-order corrections 
to this block deliver some extra factors of a and Za. In other words, a recoil 
correction with a relativistic nucleus is of the second order of the electron-to- 
nucleus mass ratio. That fact used to be 'hidden' in the standard notation. 
When one discusses the Lamb shift, the contributions used to be presented in 
the absolute units with all (m/M) factors explicitly presented. Meantime, the 
relative units are commonly used for the corrections to the hyperfine structure 
and since the Fermi energy 

E F ~ (Zaf^-mc 2 (37) 



includes a factor of m/M, a linear correction of relative order m/M is actually 
the second order contribution in the absolute units. 

Summarizing, we note that: 

• studying the gross and fine structure and the Lamb shift in a conventional 
atom such as hydrogen, one needs to take into account the hyperfine effects 
with a relatively low accuracy; 

• calculating the gross and fine structure and the Lamb shift, one does not 
need to learn much detail on theory of the second-order recoil corrections 
and to approach higher-order recoil effects; 

• accuracy of theoretical predictions suffers from an uncertainty in the deter- 
mination of the fundamental constants such as Rydberg constant Roo, the 
fine structure constant a and the nuclear mass and magnetic moment. 

Positronium is an exceptional atom since m/M = 1 and the situation strongly 
differs from the conventional atomic systems: 

• studying the gross and fine structure of positronium, one needs to take into 
account the hyperfine effects with a high accuracy; 

• the calculation of truly recoil corrections with relativistic nucleus is equally 
important for the gross and fine structure, as well as for the hyperfine struc- 
ture and the lifetime calculations. Actually, the uncertainty of the calcula- 
tions of the recoil effects determines the entire theoretical uncertainty for 
the most wanted spectroscopic characteristics in positronium. 

The equality of m/M to unity leads to another important feature of positro- 
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nium tests for the bound state QED. Since the corrections of interest are 
enhanced (m/M is not a suppressing factor any longer), the fractional accu- 
racy for successful high-precision tests is now relatively low. As a result, in 
contrast to hydrogen, the interpretation of the measurements of the IS — 2S 
interval does not crucially involve knowledge of the Rydberg constant with 
high accuracy A study of the hyperfine interval does not require a value of 
the fine structure constant a with high accuracy as it is in muonium. Since 
m/M = 1, it is unnecessary to determine a mass ratio m/M and a ratio 
of magnetic moments of the electron and nucleus, i.e., positronium 8 , in an 
additional experiment. In other words, positronium offers several high preci- 
sion tests of bound state QED without any need to determine values of any 
fundamental or phenomenological constant with high accuracy. 

I 1 1 1 1 1 1 1 1 1 1 1 1 



a i — • — i 



b i • 1 

Theory 

i i i i i i i i i i i 

203 385 203 390 203 395 

Hyperfine interval in positronium [MHz] 

Fig. 23. The hyperfine splitting in the ground state of positronium. The Yale ex- 
periment (a) was performed in 1984 [170] and the Brandeis one (6) in 1975 [171]. 

The HFS interval has been determined experimentally with the highest abso- 
lute accuracy among other positronium transitions. Present experimental data 
for the positronium hyperfine interval together with the theoretical result are 
depicted in Fig. 23. The theoretical status for the positronium hyperfine inter- 
val in the 15* state is presented in detail in Table 13. The dominant source of 
the uncertainty is the third-order corrections for positronium (AE(QED3)), 
which are known only in the logarithmic approximation. The diagrams under 
question are related to radiative-recoil and recoil corrections, and they are the 
same as those responsible for the uncertainty of muonium hyperfine interval 
(a(Za) 2 m/M and (Za) 3 m/M) [61,168,114] (see the previous section for de- 

8 One may wonder what should happen if the CPT symmetry is broken and m e /me 
and n e / are not exactly equal to unity. In principle, certain features of positronium 
might be sensitive to such a violation, however, special experiments with selective 
sensitivity to this effect would be needed. It is also clear, that such a CPT violation 
should show itself not just as a small mismatch in masses and magnetic moments, 
but as incorrectness of the basic equations (e.g., the Dirac equation). 
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tail) . The fourth-order recoil terms of theory of heavy atoms become the third 
order corrections in positronium. 



Term 




Fractional 


AE 


References 






contribution 


[MHz] 




Ep 


4 2 

a mc 


1.000 000 


204 386.6 




QED1 


a mc 


- 0.004 919 6 


-1005.5 




QED2 


fi 2 

a mc 


0.000 057 7 


11.8 


[184] 


QED3 


7 2 

a mc 


- 0.000 0061(22) 


- 1-2(6) 


[61,149,150,185] 


Total 




0.995 1321(22) 


203 391.7(6) 





Table 13 

Theory of the IS hyperfine interval in positronium. The numerical results are pre- 
sented for the frequency E/h. The calculation of the second order terms was com- 
pleted in Ref. [184], the leading logarithmic contributions were found in [61], while 
the next-to-leading logarithmic terms were achieved in [149,150,185]. The uncer- 
tainty is presented following [168,169]. 

Studies of the spectrum and decay rates of positronium provide us with a 
number of strong tests of bound state QED, some of which are among the most 
accurate. Some theoretical predictions from Table 12 can be compared with 
accurate experimental data, a review of which can be found in Ref. [166]. The 
most accurately measured spectroscopic data are related to the ground state 
HFS (see Fig. 23) and to the IS-2S interval (see Fig. 24 and Table 14). There 
are some minor discrepancies between experimental and theoretical data. 



b 




Theory 

i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 

1 233 607 200 1 233 607 220 1 233 607 240 

Positronium 1s - 2s interval [MHz] 

Fig. 24. Determination of the l 3 5i — 2 3 Si interval in positronium. The references 
here are: a is for [172] and b is for [173]. 

The experimental situation with the orthopositronium decay (see Fig. 25) 
had been not acceptable for a while but it has been recently improved [182]. 
The problem was a significant inconsistency of various experiments and long 
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Term 


A TP 




[IVirizJ 


2 2 

ol mc 


1 OQQ Con 7Q 1 

i zoo t>yu / oo. 1 


4 2 

a mc 


-8z UU5.D 


5 2 
Q TTiC 


1 /i 

-1 0U1.4 


fi 2 

a mc 


-7.1, [174] 


7 2 

a mc 


1.2(6), [175,200] 


Total 


1233 607 222.2(6) 



Table 14 

Theory of the l 3 <Si — 2 3 Si interval in positronium. The corrections AE are presented 
in the energy units, while their numerical values are given in the frequency units 
for E/h. 

standing strong discrepancy between theory and the most accurate experimen- 
tal data. The data presented in Fig. 25 include most recent vacuum results 
from Tokyo [183] and Ann Arbor [182] which are in a good agreement with 
theory (see Table 15). The original gas result from Ann Arbor (data point c 
in Fig. 25) is corrected (data point d) according to the preliminary analysis 
in Ref. [166] but that is not a final result. Further examination is in progress 
and it seems that the final uncertainty will be bigger [179]. 



□ 


Vacuum 


l — • 


— 1 h 




Gas 


l — 


□ — 1 9 


• 


Powder 


I m 






Theory 

7.025 7.035 7.045 7.055 

Decay rate of orthopositronium [u.sec-1 ] 

Fig. 25. Measurements of the annihilation decay rate of orthopositronium: a - [176], 
b - [177], c - [178], d - [166,179], e - [180], / - [181], g - [182], h - [183]. 

For the parapositronium decay theory and experiment are in perfect agree- 
ment (see Fig. 26 and Fig. 27). Most of positronium decay experiments were 
reviewed in Refs. [166,167]. 

These two papers also review experiments on the fine structure in positronium 
performed at 2 3 5*i — 2P intervals which were less accurate than experiments at 
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Contribution 


Decay rate 


Decay rate 




of orthopositronium 


of parapositronium 




r In 

[^sec J 


r In 

[/usee J 


r(°) 


7.211 17 


8 032.50 


QEDl 


-0.172 30 


-47.25 


QED2 


0.001 11(1), [186] 


4.43(1), [187] 


QED3 


-0.000 01(2), [61,188] 


-0.08(4), [61,188] 


Total 


7.039 96(2) 


7989.62(4) 



Table 15 

Theory of the annihilation decay rate of ortho- and parapositronium (the 15 state). 
The leading contributions are defined above in Table 12. The decay rate of or- 
tho/parapositronium into five/four photons is included into corresponding QED2 
terms. 





i ' 

Ann Arbor, 
vacuum 




— • — I 


Theory 

i 



7 980 7 990 8 000 

decay rate 
of parapositronium [fisec" 1 ] 



Fig. 26. Determination of the parapositronium decay rate. The experiment was 
performed in Ann Arbor in 1994 [189] . 



d h 



a h 



H e 



1 c 

b 



1.0 1.5 2.0 

Parapositronium branching 4y/2y [ppm] 



Fig. 27. Annihilation of parapositronium: branching Br(47/27) measured around 
the world (a - [190], b - [191], c - [192], d - [193], e - [194]). 
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IS HFS and IS — 2S intervals. Data on the fine structure are not as accurate 
as the results related to the IS" hyperfine splitting and the IS — 2S interval 
but the progress is possible [166]. 




8620 8630 8640 11 170 ' 11 180 ' 11 190 13 000 ' 13 010 ' 13 020 18 500 18 510 18 520 

Positronium 2S-2P fine structure [MHz] 

Fig. 28. The fine structure in positronium at n = 2: the energy levels and studied 
transitions. The 2 3 S\ — 2 l P\ transition is forbidden for the free positronium atom 
and it was measured with an external magnetic field applied. The line width of the 
2 3 Si state is of 1 MHz and that of all 2P states is of 50 MHz and thus it cannot 
be seen in the scale of the top picture, while the width of parapositronium singlet 
(2 1 Sq) being of 1.3 GHz is clearly recognized. The references to the measurements 
are: a - [195], b - [196], c - [197], d - [198] and e - [199]. 

To conclude the section, let us discuss briefly the theory of the positronium 
decay rates and energy levels and in particular the uncertainty of their cal- 
culations. We collected all theoretical predictions in Table 12. The results 
were published and presented in different compilations. What we would like 
to emphasize here is our estimation of uncertainty . 

The uncertainty of any perturbative theoretical calculation is determined by a 
plausible estimation of unknown higher-order corrections which are expected 
to have large coefficients. There is a number of corrections enhanced by a big 
double logarithmic factor In 2 a ~ 24 [61] and the higher-order terms should be 
studied to understand better the accuracy of theory. For the spectroscopy, the 
crucial order is a 7 mc 2 , while in the case of the decay the theoretical problems 
are related to the corrections of the relative order a 3 . 

For most of them not only the leading logarithmic corrections (e.g., in the case 
of the spectrum that is a 7 mc 2 In 2 a) are known, but also the next-to-leading 
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LillivllU WI1 
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\~\ 7/11 

[174J 


[1 1 0,ZUUJ 


unknown 


15 HFS 


[184] 


[61] 


[149,150,185] 


r(p-Ps) 


[187] 


[61] 


[188] 


r(o-Ps) 


[186] 


[61] 


[188] 



Table 16 

References to recent progress in the positronium theory. The contributions to the 
spectrum are classified by the absolute energy units (such as a 6 mc 2 ), while those 
to the decay are presented in units of the leading contribution T (see Table 12). 

term (a 7 mc 2 In a). We should emphasis, however, that such a knowledge can- 
not reduce the uncertainty. The leading term originates (the most probably, 
but not always) from a single source and its magnitude is a natural charac- 
teristic of the correction, while the next-to-leading term used to be a result of 
a summation of or a cancellation between different contributions and can be 
sometimes small. But that smallness is misleading and the constant following 
the single logarithm is not necessary small. Our estimation of the uncertainty 
[168,169] is based on a value of the double logarithmic term in any case. 



10 The g factor of an electron and a nucleus in a light hydrogen-like 
atom 

10.1 General consideration 

Not only a spectrum of simple atoms can be studied with a high accuracy. 
Among other quantities accessible to high precision measurements are the 
magnetic moments of bound subatomic particles. The only way to determine 
a value of a magnetic moment precisely (i.e., with uncertainty at the level 
of 0.01 ppm or below) is to measure a certain frequency proportional to the 
magnetic moment and applied magnetic field. Since the magnetic field can be 
neither measured nor calculated with uncertainty essentially below 0.1 ppm, 
it is possible to measure accurately only certain ratios of the magnetic mo- 
ments. Another option is to compare a spin precession frequency related to the 
magnetic moment with the ion cyclotron frequency. Two-body atomic systems 
provide us with an opportunity to study the magnetic moments of their bound 
composites and in certain experiments to compare their magnetic moments to 
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each other. 



The interaction of an atom with a weak homogenous magnetic field B can be 
expressed in terms of an effective Hamiltonian. For an s state in the two-body 
atom such a Hamiltonian for spin variables has the form (see, e.g., [9]) 

# mag n = 9'e 77- (s • B) - g' N — T (I • B) + J HFS (s • I) , (38) 



where s stands for the spin of an electron, I is for the nuclear spin, g' e t N \ is 
for the g factor of a bound electron (nucleus) in the two-body atom and Jhfs 
is the hyperfine constant 9 . We remind here, that m and M are the masses of 
the electron and the nucleus, respectively. 

We present the g factor of a bound electron in the form [60,202,203] 

g' = 2 • (l + a + b) , (39) 



where a is the anomalous magnetic moment of a free electron, while b stands 
for the binding correction. As it has been known for a while [9,204], even in 
absence of QED effects the g factor of a bound electron is not equal to the 
free Dirac value (namely, g^l = 2). The leading radiative and recoil effects 
were later examined in [205,206]. 

Three basic situations with the two-body atomic systems at presence of a 
homogeneous magnetic field are of practical interest: 

• a 'conventional' low-Z atom with a non-zero nuclear spin (hydrogen, deu- 
terium, tritium, helium-3 ion); 

• muonium, a pure leptonic unstable atom with a relatively light nucleus 
(m e /m M ~ 1/207 > m e /m p ~ 1/1836); 

• a medium- Z hydrogen-like ion with a spinless nucleus (the ions of carbon- 12 
and oxygen-16) with a possible expansion of experimental activity down to 
helium-4 and beryllium- 10 ions and up to calcium-40 ion. 

The bound g factors in light atoms are now known up to the fourth-order 

9 The notation for the g factors is different in particle and nuclear physics. For ions 
one needs to apply a different notation in which the magnetic Hamiltonian reads 

/W = ^£(s-B)-<4 — (I-B) + J HFS (s-I), 

where g' N differs indeed from that in Eq. (38). To avoid confusion, in this paper g is 
used only for the case of particles (electrons, muons, protons etc.), while if we need 
the g factor for a nucleus, we introduce additional notation (see e.g., (cf. 34)). 
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B 



Fig. 29. Sublevels E maLga (F, F z ) of the hyperfine structure in the ns state of hydrogen 
or another hydrogen-like atom with nuclear spin 1/2 at a homogenous magnetic field 
B. The related equation was first derived by Breit and Rabi (see [9]). 

corrections [207] including the terms of the order a 4 , a 3 m/M and a 2 (m/M) 2 
and thus the relative uncertainty is substantially better than a part in 10 8 . 
All these corrections are of kinematic origin and there is no need for a higher 
accuracy and thus for any higher order corrections for a moment. 

For the medium- Z ions some corrections of higher orders are important for the 
comparison with experimental data and for applications to the determination 
of the fundamental constants. These corrections involve essential bound state 
QED effects. The 'kinematic origin' reads that they are a result of a 'simple' 
quantum mechanical problem of two point-like particles bound by a Coulomb 
interaction with their actual values of the masses and magnetic moments. 
'Free QED' corrections are a part of kinematic effects since we consider the 
actual values of the magnetic moments which include the anomalous magnetic 
moment. 'Bound state QED' corrections are related to somewhat beyond this 
approximation — to non-Coulomb additions to the interactions (as, e.g., the 
Uehling potential), structure effects (as, e.g., the electron self-energy). Such 
higher order corrections were studied and in particular the linear correction 
in a has been known exactly in all orders of Za as well as some higher order 
recoil contributions. 



10.2 Hydrogen and its isotopes: the isotopic effects for the g factor of a bound 
electron 

Few experiments were performed with hydrogen and its isotopes, as well as 
with the helium-3 ion. Dependence of energy sublevels on the magnetic field 
is presented in Fig. 29 for the ns state in a two-body atom with nuclear spin 
1/2 (hydrogen, tritium, helium-3 ion, muonium). 
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Value Experiment Reference Theory 



5(H) /g (e) 


1 


- 17.709(13) • 10~ 6 


[208,209] 


1 


- 17.694 • 10~ 6 


<7(D)/<?(H) 


1 


- 7.22(3) • 10~ 6 


[210] 


1 


- 7.24 • 10" 6 


<7(T)/<?(H) 


1 


- 10.7(15) • 10- 6 


[211] 


1 


-9.7- 10~ 6 


ff ( 4 He+)/ 50 (e) 


1 


- 70.88(30) • 10" 6 


[208,209,212,213,214,215] 


1 


- 70.91 • 10~ 6 



Table 17 



The bound electron g factor in light hydrogen-like atoms: a comparison of theory 
and experiment. Here, go{e) stands for the magnetic moment of a free electron and 
it contains the anomalous magnetic moment a = a/2ir + ... 

The hyperfine sublevels were successfully studied by a few different ways. 

• For example, one is a study of a few splittings related to the same atom. 

• Other experiments on the hydrogen-like atoms dealt with the g factors of 
an electron bound in different atoms studying a kind of isotopic dependence 
in the electron g factor. 

• One more approach was realized for muonium: the field-induced energy 
intervals were measured at the magnetic field calibrated against the proton 
spin precession (see Sect. 10.4). 

The isotopic experiments were performed on hydrogen, deuterium and tritium 
measured at the same applied magnetic field. The results are summarized in 
Table 17. The helium- ion result there involved a chain of several intermediate 
comparisons including experiments with neutral rubidium. All corrections to 
the electron g factor up to the fourth order are known, with an exception. 
There are two contradicting results for the (Za) 2 (m/M) 2 term which can be 
potentially nuclear-spin-dependent. One group [216] performed a semi-classical 
evaluation and reported on a nuclear-spin-independent result, while the others 
[217] claimed a difference between the results for the nuclear spin 1 = 1/2 
and 1=1. For the nuclear spin 1/2 both reproduced the well known result 
[206] . Hopefully, for deuterium their discrepancy is somewhat smaller than the 
experimental uncertainty, while the other two-body atoms of interest involve 
the nuclear spin /= 1/2 for which the result in this order is well established. 

Setting the nuclear spin equal to 1/2, the equations for the bound g factor in 
two-body atom read [207,218] 
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97 
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where the anomalous magnetic moment of the nucleus is defined as 



where, we remind, M is the nuclear mass, / is its spin, while /xjv m the nuclear 
magneton. Thus a p ~ 1.793..., ~ —0.143..., a t ~ 7.916... and ah ~ 
— 1.184. . . (the latter value is for a helion, the nucleus of the 3 He). We note 
that appearance of the factor of (1+ajy) in the denominator is rather artificial: 
a more natural (from the theoretical point of view) expression is of the form 

g' N = 2 x (1 • g 1 + a N ■ g a ) , 



however, due to applications is it more convenient to present the result as 
above (i.e., as a multiplicative correction to gx, a free value of the g factor). 

We remark also that some of the listed values of the nuclear anomalous mag- 
netic moments are of large numerical values (as, e.g., a t ). That is not an 
exception. The triton magnetic moment is essentially a proton magnetic mo- 
ment (the total spin of two neutrons is zero) somehow perturbed. However, the 
magneton related to the triton is approximately a third part of the proton's 
since the charge of the proton and triton is the same, but the triton mass is 
threefold higher. 



10.3 Hydrogen and deuterium: determination of the proton and deuteron 
magnetic moments 



Studies of a few HFS intervals in the same atom in the presence of the magnetic 
field were performed for hydrogen [219] and deuterium [220] and their results 
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2m 




(42) 
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were found to be 



A* e (H) _ ff e (H) hb 
/i p (H) g p (R) jj, N 



658.210 705 8(66) 



(43) 



and 



A*e(D) 



2143.923 565(23) . 



(44) 



The bound state QED contributions are small and the quoted experiments 
were not designed to test them. The purpose was different, namely, to de- 
termine the proton magnetic moment which is widely used as a probe mag- 
netic moment for other QED tests (such as the muonium hyperfine structure 
[131,135], the anomalous magnetic moment of a muon [137]) or for other mea- 
surements (such as of the deuteron magnetic moment [129] etc.). In fact, the 
hydrogen experiment [219] with the theoretical evaluation performed in [218] 
has delivered the most accurate value of a free ratio of the electron-to-proton 
magnetic moments 



and the proton g factor. The result of the evaluation in [218] is somewhat 
shifted from the CODATA value [21] due to recently obtained higher-order 
corrections [221,222,207] and improvement in determination of the electron- 
to-proton mass ratio discussed below (see Sect. 11). 

In the case of deuteron there are competitive NMR experiments on HD spec- 
troscopy. The experimental results published about twenty years ago in [223,224,225,226] 
were obtained in NMR high-pressure experiments and they are not free of pos- 
sible systematic effects. We start a program of additional experiments [129] 
and we hope to study those systematic effects in more detail and achieve 
reliable data. The present situation is summarized in Fig. 30. 



Here and further in this section we ignore the direction of the spin and the mag- 
netic moment and thus the sign of some g factors and ratios of magnetic moments. 
To simplify notations, we also drop the prime sign when the bound system is ex- 
plicitly specified as, e.g., in the case of g(H). 



He 



658.210 685 9(66) 



(45) 
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Fig. 30. The deuteron-to-proton ratio of the magnetic moments: maser result a is 
taken from [219,220,21], the other results are from NMR spectroscopy of hydrogen 
deuteride HD. Result b is related to an experiment with an admixture of helium-3 
at pressure of 60-80 atm [225], while results c and d correspond to measurements on 
pure HD gas at pressure of 130 atm with/without rotation of the sample [226]. The 
high-pressure results here have an uncertainty corrected for an inaccuracy in the 
calculation of screening effects (see [129] for detail). The result e is a recent result 
obtained at a medium pressure (10 atm) in order to check effects of the pressure 
dependence [129]. The systematic effects due to the pressure shift are under study 
and have not been included into the uncertainty budget. 

10. 4 The g factors of an electron and a muon in muonium and the determi- 
nation of the muon mass 



A difference between theory for muonium and for 'conventional' atoms is due 
to a higher importance of recoil effects (since the nucleus, namely a muon, is 
approximately tenfold lighter than a proton) and a possibility to expand the 
muon magnetic moment in a. After expanding a M , the result for the bound 
muon g factor reads [207,114] 
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_ „(Z,)(1 + Z) / my g^gdO m _ 97 3 1 (47) 

2 \M7 12tt M 108 V ; J ' V ; 

where the g factor of a free lepton includes the anomalous magnetic moment 
^2 = 2.(1 + 0^). 

Studies of the Breit-Rabi levels (see Fig. 29) in muonium delivered the best 
datum on the electron-to-muon mass ratio needed for a determination of the 
fine structure constant a from the muonium hyperfine structure (see Sect. 8). 
Equation (38) has been applied [131,135] to determine the muon magnetic 
moment and the muon mass by measuring the splitting of sublevels in the 
hyperfine structure of the Is state in muonium in a homogenous magnetic field. 
Since the magnetic field was calibrated via the spin precession of a proton, 
the muon magnetic moment was measured in units of the proton magnetic 
moment, and the muon-to-electron mass ratio was derived as 

m v _ /V ih_ \ ug\ 



The results on the muon mass extracted from the Breit-Rabi formula are 
among the most accurate (see Fig. 31). Still there is a more precise value 
available derived from the muonium hyperfine structure after a comparison 
of the experimental result with theoretical calculations. However, the latter 
is of reduced interest, since the most important application of the precise 
value of the muon-to-election mass is to use it as an input for calculations 
of the muonium hyperfine structure while testing QED or determining the 
fine structure constants a. The adjusted CODATA result [21] in Fig. 31 was 
extracted from the muonium hyperfine structure studies and in addition used 
a somewhat over-optimistic estimation of the theoretical uncertainty (see for 
detail [133]). 



11 The g factor of a bound electron in a hydrogen-like ion with a 
spinless nucleus and a determination of m e /m p 



A study of the g factor of a bound electron in a medium-Z hydrogen-like 
ion with a spinless nucleus offers a comparison of the electron spin precession 
frequency affected by QED effects and the ion cyclotron frequency. That can 
be used to test QED theory as well as to determine the electron-to-proton 
mass ratio. An important feature of such a study is that in contrast to most of 
the other tests it is possible to do both without any interference between these 
two tasks. That is possible because of an opportunity to perform experiments 
with different hydrogen-like ions. 
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Fig. 31. The muon-to-electron mass ratio. The most accurate result obtained from 
comparison of the measured hyperfine interval in muonium [131] to a theoretical cal- 
culation [133] performed with oT g \ = 137.035 998 76(52) [134]. The results derived 
from study of the Breit-Rabi sublevels are related to two experiments performed at 
LAMPF in 1982 [135] and 1999 [131]. The others are taken from the measurement 
of the Is — 2s interval in muonium [130], muon-spin-precession-resonance (MSR) 
study of a muon in bromine [136] and from the CODATA adjustment [21]. 

Since the anomalous magnetic moment of a free electron is known with high 
accuracy [227,21] 



a e = 0.001159 652 187(4) 



(49) 



to calculate the bound g factor denned in Eq. (39), one has to obtain a bound 
correction b. A summary of the calculation of the bound corrections is pre- 
sented in Table 18. We mainly follow here our consideration in [207] 11 . In 
particular, the uncertainty of unknown two-loop contributions is taken from 
[230] , where the theoretical uncertainty was reasonably estimated for the first 
time. The results of the one-loop self-energy are achieved in different ways for 
different atoms. For lighter elements (helium, beryllium) it is based on fitting 
[207] data of [231], while for heavier ions we use the results of [232]. The other 
results are taken from [230] (for the one-loop vacuum polarization; see also 
[233]), [203] (for the nuclear correction and the electric part of the light-by- 
light scattering (Wichmann-Kroll) contribution; see also [233]), [234] (for the 
magnetic part of the light-by-light scattering contribution) and [235] (for the 
recoil effects). 

11 Recently a certain progress in the calculation of the one-loop and two-loop self- 
energy contribution [228] as well as of the the magnetic part of the light-by-light 
scattering contribution [229] was achieved. Since presently the experimental uncer- 
tainty dominates, this improvement of theory has no consequences for a moment. 
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Ion g 

4 He+ 2.002177 406 7(1) 

10 Be 3+ 2.001751574 7(2) 

12 C 5+ 2.001041590 0(4) 

16 7+ 2.000 047020 3(8) 

18 7+ 2.000 0470211(8) 

Table 18 

The bound electron g factor in low-Z hydrogen-like ions with spinless nucleus. 

Before comparing the theory to experimental data, let us shortly describe 
certain details of the experiment crucial for its interpretation. To determine 
a quantity like the g factor, one needs to measure a certain field-induced 
frequency (such as the spin precession frequency) at some known magnetic field 
B. It is clear that there is no way to determine the strength of the magnetic 
field with a high accuracy directly. The conventional way is to measure two 
frequencies and to compare them. The frequencies measured in the GSI-Mainz 
experiment [236] are the ion cyclotron frequency 

= [ ~M^ B • < 50) 

where Mj is the ion mass, and the Larmor spin precession frequency for a 
hydrogen-like ion with a spinless nucleus 

«L = 9b ^-B, (51) 
where gt, is the bound electron g factor. 

Combining them, one can obtain a result for the g factor of a bound electron 

or for the electron-to-ion mass ratio 

^ = ^. (53) 
Mi Z -1 2 u L K ' 



To measure the frequencies, one has to apply a certain magnetic filed and, in 
principle, it should be as homogenous as possible. However, the measurement 
of a frequency consists of two important parts: an induction of a certain mag- 
netic transition by an accurately calibrated perturbation and a detection of 
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the transition. The homogeneity of the field, which is a crucial requirement of 
the induction, is in a contradiction with the detection scheme, which requires 
a non-homogeneity of the magnetic field. In particular, a certain z-dependence 
of the field was applied providing ion oscillations in a vertical direction which 
were detected. Since the oscillation frequency depends on the magnetic state 
that allows to identify the state and thus to detect the transition. The break- 
through in a study of the bound g factor [236] was a result of a space separation 
of the transition and detection areas, which gives an opportunity to support 
a very homogenous field for the transitions and afterwards to measure the 
response function in a separate area with z gradients. 

Today the most accurate value of m e /Mj (without using data from the experi- 
ments on the bound g factor) is based on a measurement of m e /m p realized in 
a Penning trap [237] with a fractional uncertainty of 2 ppm. The accuracy of 
measurements of u c and u>l as well as of the calculation of g^ (as shown in [202]) 
is substantially better. That means that it is preferable to apply Eq. (53) to 
determine the electron-to-ion mass ratio [238]. Applying the theoretical value 
for the g factor of the bound electron and using the experimental results for 
uj c and ujl in hydrogen-like carbon [236] and some auxiliary data related to 
the proton and ion masses from [21], we arrive at the following values 

—^ = 1836.152 6731(10) (54) 
m e 

and 

m e = 0.000 548 579 909 29(31) u , (55) 

which slightly differ from those in [238]. The present status of the determina- 
tion of the electron-to-proton mass ratio is summarized in Fig. 32. 

In [202] it was also suggested in addition to the determination of the electron 
mass to check theory by comparing the g factor for two different ions. In such 
a case the uncertainty related to m e / Mi in Eq. (52) vanishes. The theoretical 
uncertainty is chiefly related to uncalculated higher-order terms. Since various 
higher-order corrections are proportional to Z 4 or so, even a small change in 
Z, such as, e.g., a comparison of the carbon result [Z = 6) to the oxygen one 
(Z = 8) is a quite sensitive probe for the value of the higher-order terms. 

Thus, combining the experimental results for carbon [236] and oxygen [239] 
ions 

#( 12 C 5+ )M 16 7+ ) = 1.000497273 1(15) , (56) 
one can verify the reliability of the estimation of the theoretical accuracy. The 
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Fig. 32. The proton-to-electron mass ratio. The theory of the bound g factor is taken 
from Table 18, while the experimental data on the g factor in carbon and oxygen 
are from [236,239]. The Penning trap result from University of Washington is from 
[237]. 

experiment appears to be in perfect agreement with theory 

#( 12 C 5+ )M 16 7+ ) = 1.000 497273 3(9) . (57) 

Unfortunately, the authors of the experiments did not present this value di- 
rectly [239] . They presented results for the g factors of both ions separately. 
In both cases the experimental uncertainties contain systematic and statisti- 
cal contributions. The dominant contribution is due to the electron mass and 
it should vanish in the ratio. It is also not good to determine m e /Mi from 
one experiment (let us say from carbon's) and use as an input for the other 
(oxygen's). In such a way additionally to the experimental uncertainties the 
experimental g factor for the oxygen would contain a theoretical uncertainty 
for the carbon g factor. 

That is why it is so important to speak in terms of the ratio of the g factors. 
The experimental ratio is to be determined by experimentalists who perfectly 
know the correlations between the systematic sources of two measurements, 
while the theoretical ratio is to be determined by theorists and that allows 
to take into account correlations in the estimations of unknown higher-order 
terms. 

The estimation in Eq. (56) is based on our understanding of the papers 
[236,239] and on a private communication [240]. The agreement between the- 
ory and experiments means that we have a reasonable estimate of the uncal- 
culated higher-order terms. Note, however, that for metrological applications 
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it is preferable to study lower Z ions (hydrogen-like helium-4 and beryllium- 
10) to eliminate these higher-order terms if the experimental accuracy will 
increase. 

As mentioned above, the leading radiative corrections are of kinematic origin 
and in particular the contributions of order a(Za) 2 and a 2 (Za) 2 are related 
to a relativistic bound-state Coulomb problem for an electron and muon with 
non-zero anomalous magnetic moments, which are determined by free QED. 
The leading bound-state QED effects are of order a(Za) 4 (the one-loop con- 
tributions) and a 2 (Za) 4 (the two- loop contributions). The former are known 
with sufficient accuracy, while lack of accurate results for the latter determines 
the uncertainty of the theory. That is compatible with a test of the Lamb shift 
theory at the level of a percent (for the hydrogen-like ions of carbon and oxy- 
gen). That is indeed not accurate enough for any interesting bound-state QED 
test. And in fact that is a big advantage of the experiments on the bound elec- 
tron g factor for determination of m e /m p . In a typical situation (the Lamb 
shift and the Rydberg constant; the muonium hyperfine structure and the fine 
structure constant; the helium fine structure and a) we have to try to verify 
QED calculations and to determine a certain fundamental constant within the 
same experiment. In the case of the determination of the electron-to-proton 
mass ratio the QED theory is relatively simple. We still need to prove QED 
calculations, but we have not yet approached the crucial problem of the Lamb 
shift theory, namely, a problem of higher-order two-loop corrections. However, 
with sufficient theoretical progress this problem can be studied with present 
experimental accuracy if the experiment will be turned to higher Z, in par- 
ticular, to Z = 20 — 30. The higher-order bound-state QED effects should 
contribute at a detectable level and may be successfully studied. Comparing 
these higher Z results with Z = 6, 8 and lower Z one can cancel an uncertainty 
due to the determination of the electron-to-proton mass ratio. Thus, the ex- 
perimental verification of bound-state QED theory and the determination of 
m e /m p can be clearly separated. 



12 A determination of the fine structure constant a by means of 
QED and atomic physics 

Precision measurements accompanied by accurate theoretical calculations strongly 
interfere with high precision determinations of certain fundamental constants. 
Those are needed in order to interpret the theoretical expressions in terms 
of measurable quantities. The fine structure constant a plays a basic role in 
QED tests. In atomic and particle physics there are several ways to deter- 
mine its value. The results are summarized in Fig. 33. One method based on 
the muonium hyperfine interval was briefly discussed in Sect. 8. A value of 
the fine structure constant can also be extracted from the neutral-helium fine 



72 



structure [241,242,243] and from a comparison of theory [134] and experiment 
[227] for the anomalous magnetic moment of an electron (a g _2)- The latter 
value has been the most accurate one for a while and there was a long search 
for another competitive value. Since the value of a s _ 2 is based on one-group 
experiment [227] and one-group theory [244] there has been a long stand- 
ing concern on reliability of this value and a search for its confirmation, in 
particular, via an independent competitive determination of the fine structure 
constant. At present, the second best value (ctcs) m the list of the most precise 
results for the fine structure constant comes from recoil spectroscopy [245]. It 
is only three times less accurate than a 9 _2- We note, however, that the result 
has not yet been published in a refereed journal, but only in several confer- 
ence proceedings [245]. Since results in proceedings are less reliable than those 
published in the journals, we look forward to the eventual publication. 



helium fs (v 2 i 
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Cs recoil 



CODATA, 1998 



137.035980 137.036000 137.036020 

QED values of the inverse fine structure constant a 1 



Fig. 33. The value of the fine structure constant a as determined by means of atomic 
physics and QED. 

We would like to consider briefly the use and the significance of the recoil 
result for the determination of the fine structure constant. Absorbing and re- 
emitting a photon, an atom, if it is initially in rest, gains some kinetic energy 
which can be determined through a measurement of the shift of the emitted 
frequency in respect to the absorbed one. The recoil shift depends on relative 
directions of two photons 



Af = -C{6) 



hf 2 
2M a c 2 



(5c 



where M a is the mass of the atom, and under conditions of the experiment 
[245] C = 2. 



That is a measurement of the frequency Af with high accuracy that was the 
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aim of the photon recoil experiment [245] . Combining the absorbed frequency 
and the shift Af, it is possible to determine a value of atomic mass in frequency 
units, i.e., a value of M a c 2 /h. In experiment [245] the atom was cesium and 
the transition was the Dl line. The atomic mass of cesium is known very 
well in atomic units (or in units of the proton mass) [246] and the Dl line 
was accurately measured in [247]. A value of M a c 2 /h may be compared to 
the Rydberg constant Roo = a 2 m e c/2h and if the electron mass in proper 
units is known, that eventually delivers us a value of a. A determination 
of the electron mass is more complicated than of any atomic mass because 
of a gap of three orders of magnitude between masses of an electron and a 
proton. At present, the biggest uncertainty of the recoil photon value of ac s 
comes from the experiment [245], however, accuracy in the determination of 
the electron mass is still significant since that is the second important source 
of the uncertainty. 

The success of the ac s determination was sometimes ascribed to the fact that 
Ug-2 is a QED value, derived with the help of QED theory of the anoma- 
lous magnetic moment of an electron, while the photon recoil result is free of 
QED. We would like to emphasize that the situation is not so simple and the 
importance of the involvement of QED is rather exaggerated. What is more 
important is that the uncertainty of u g -2 originates from understanding of 
the electron behavior in the Penning trap and it dominates over any QED 
uncertainty. For this reason, the value of ac s , based on m p /m e from another 
Penning trap experiment [237] and actually obtained by the same group, which 
measured the value of the anomalous magnetic moment of electron [227] , can 
actually be correlated with a g _ 2 . The result 

acl = 137.036 0002 8(10) (59) 



presented in Fig. 33 is obtained using m p /m e from Eq. (54). The value of the 
proton-to-electron mass ratio found this way is free of the problems with an 
electron in the Penning trap, but a certain QED theory is involved. However, 
one has to realize that the QED uncertainty for the g factor of a bound 
electron and that for the anomalous magnetic moment of a free electron are 
completely different. The bound state QED theory deals with relatively simple 
Feynman diagrams but in Coulomb field and in particular to improve theory 
of the bound g factor, we need better understanding of Coulomb effects for 
'simple' two-loop QED diagrams. In contrast, for the free electron no Coulomb 
field is involved, but the uncertainty of the computation arises because of the 
accuracy in a numerical calculation of a big number of complicated four-loop 
diagrams. There is no correlation between these two kinds of calculations. The 
words 'simple' and 'complicated' are related to the diagrams, but not to the 
simplicity and difficulties of the calculations. Both kinds of calculations are 
quite difficult but in a different way. The crucial free diagrams are complicated 
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diagrams built of simple blocks. The key bound QED contributions correspond 
to a simple diagram constructed of complicated blocks. 

The two best values, cx g -2 and ac s , agree with the muonium result but contra- 
dict to the helium value (see Fig. 33). The latter is presented in Fig. 33 rather 
to show a potential accuracy of the application of the helium fine structure 
method. The helium fine structure of the ls2p states involves few transitions 
and in particular so-called 'small interval' (z^oi) and 'big interval' (^12) • A com- 
parison of theory against experiment is presented in Fig. 34. One can match 
a certain scatter of the experimental data and a bad agreement of theory 
and experiment. In principle, one could speculate on a special value of the 
fine structure constant which will set an agreement of theory and experiment. 
However, the contradiction for z/ i and u\ 2 is at the same level in absolute 
units (from 5 to 20 kHz) and differs in fractional units since the 'big interval' 
is more than tenfold bigger than the 'small interval'. The contradiction for 
the 'small interval' is too big to be explained by any acceptable shift in the 
value of a. The situation is quite unclear but, still, the recent improvement of 
theoretical accuracy remains quite promising [241,242]. 

There is a number of other values of the fine structure constant (see [21,24] 
for detail). Some of them (based on a measurement of the gyromagnetic ratio 
of a proton and a helion and on a determination of h/m n ) involve some QED 
theory and experiments, however, the uncertainty of all of them is completely 
determined by metrological problems due to electric standards and a produc- 
tion of a perfect silicon crystal for the Avogadro project. The other electrical 
determinations, based on a measurement of the quantum Hall resistance in 
the SI units, needs no QED. Since the accuracy of those results is substan- 
tially below that of a 9 _ 2 and since their uncertainty is related to the electric, 
magnetic and material standards we do not consider them here. An updated 
review of QED and non-QED determination of the fine structure constant can 
be found in [23]. 



13 Higher-order logarithmic corrections and the uncertainty of the 
bound-state QED calculations 

13.1 Hyperfine structure in light hydrogen-like- atoms 

To summarize QED tests related to the hyperfine structure, we present in 
Table 19 the data related to the HFS interval in the Is state of positronium 
and muonium and to the specific difference D21 in hydrogen, deuterium and 
helium-3 ion. As we noted in Sects. 7, 8 and 9, their theoretical uncertainty 
is related to recoil effects of relative order a 3 (m/M). In the case of helium 
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Fig. 34. The fine structure in neutral helium, uqi stands for ls2p 3 Po — 3 -Pi an d ^12 
for ls2p 3 Pi — 3 P2- The theory is presented according to [241,242], the experimental 
data are taken a from [243], b from [248], c from [249], d from [250] and e from [251]. 
Less accurate data are not included. 



D 2 i the higher-order one-loop and two-loop effects are also important, as well 
as higher-order nuclear corrections. The theoretical accuracy supersedes the 
experimental one (except of the case of the -D21 HFS interval in the helium- 
3 ion, for which the theoretical uncertainty was recently increased after its 
reconsideration [115,116]). The theory agrees with experiment very well. A 
certain expection is the positronium case where the agreement is not perfect. 
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Atom 


Experiment 


Theory 


A/a 


a/E F 




[kHz] 


[kHz] 




[ppm] 


Hydrogen, D21 


49.13(13), [106] 


48.953(5) 


1.3 


0.09 


Hydrogen, D21 


48.53(23), [107] 




-1.8 


0.16 


Hydrogen, D21 


49.13(40), [108] 




0.4 


0.28 


Deuterium, D21 


11.280(56), [109] 


11.312 5(10) 


-0.6 


0.17 


Deuterium, D21 


11.16(16), [110] 




-1.0 


0.49 


3 He + ion, D 2 i 


-1189.979(71), [111] 


-1 190.08(15) 


0.6 


0.02 


3 He+, D 21 


-1190.1(16), [112] 




-0.01 


0.19 


Muonium, Is 


4463 302.78(5) 


4463 302.88(55) 


-0.18 


0.11 


Positronium, Is 


203 389100(740) 


203 391700(500) 


-2.9 


4.4 


Positronium, Is 


203 397 500(1600) 




-2.5 


8.2 



Table 19 

A comparison of experiment and theory of the hyperfine structure in light hydrogen- 
like atoms. The numerical results are presented for the frequency E/h. In the D21 
case the reference is given only for the 2s hyperfine interval, since for all the atoms 
of interest the Is HFS interval was measured more accurately (see table 4 for more 
detail) . 

13.2 Crucial higher-order corrections 

The precision physics of light simple atoms provides us with an opportunity 
to check higher-order effects of perturbation theory. The highest-order terms 
important for the comparison of theory and experiment are collected in Ta- 
ble 20. The results for the energy levels and decay rates of interest. There are 
also a few dimensionless quantities of interest, which are mainly the g factors. 
The uncertainty of the g factor of the bound electron in carbon and oxygen 
is related to a 2 (Za) 5 mc 2 corrections in the energy units, while for calcium 
the crucial order is a 2 (Za) 6 mc 2 if the experiment will reach the same level of 
accuracy. 

While some of the crucial corrections presented in Table 20 are completely 
known, others are not. Many of the listed corrections and in particular aiZafvn 2 
and (Za) 7 m 2 /M 3 for the hyperfine structure in muonium and helium ion, 
a 2 (Za) e m for the Lamb shift in hydrogen and helium ion, a 7 m for positro- 
nium have been known in a so-called logarithmic approximation. In other 
words, only the terms with the highest power of a 'large' logarithm (e.g., 
\n{l/Za) ~ ln(M/m) ~ 5 in muonium) have been taken into account. We 
note that the 'large' corrections should have an essentially non-relativistic ori- 
gin. Truly relativistic corrections used to have a few factors Za/n instead of 
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Value Order 

[mc 2 ] 

Hydrogen, deuterium (gross structure) a(Za) 7 , a 2 (Za) e 

Hydrogen, deuterium (fine structure) a(Za) 7 , a 2 (Za) e 

Hydrogen, deuterium (Lamb shift) a(Za) 7 , a 2 (Za) 6 

3 He+ ion (2s HFS) a(Za) 7 m/M, a(Zafm 2 /M 2 , 

a 2 (Zafm/M, {Za) 7 m 2 /M 2 
4 He+ ion (Lamb shift) a(Za) 7 , a 2 (Zaf 

Muonium (Is HFS) {Za) 7 m 2 /M 2 , a{Zafm 2 /M 2 , 

a(Za) 7 m/M 

Positronium (Is HFS) a 7 

Positronium (gross structure) a 7 

Positronium (fine structure) a 7 

Para-positronium (decay rate) a 7 

Ortho-positronium (decay rate) a 8 

Para-positronium (47 branching) a 8 

Ortho-positronium (57 branching) a 8 
Tabie~20 

A comparison of QED theory and experiment: the crucial orders of magnitude for 
the energy levels and decay rates in units of mc 2 (see [60] for more detail). 

Za, which is a characteristic non-relativistic factor for the Coulomb problem, 
and thus they are numerically suppressed. First calculations of logarithmi- 
cally enhanced higher-order terms were performed for higher-order corrections 
in muonium HFS, hydrogen Lamb shift, positronium energy level and decay 
rates in our paper [61] and were successfully developed by various authors 
in [148,34,252,175,200,124,113]. By now even some non-leading logarithmic 
terms have been evaluated by several groups [149,150,185,67,68]. 

Such a value of the parameter as ~ 1/5 is sufficient to make a good first 
approximation. We used to estimate accuracy of the leading logarithmic ap- 
proximation by 50%. However, we cannot consider such an estimation as an 
accurate and reliable result. If the correction happens to have a value signifi- 
cant in comparison with the overall uncertainty, one has to improve the result. 
A use of the leading logarithmic terms is fruitful since they used to have some 
'natural' values because they are always (or almost always) state-independent 
and not a result of big cancellations. The latter happens seldom between differ- 
ent gauge-invariant sets (e.g., in the aEp term for the positronium hyperfine 
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splitting) and is easy to recognize. It is clear what to do in the case of such 
a cancellation; it is necessary to estimate the non-leading terms separately by 
a half-value of the related leading terms and then to sum these uncertainty 
contributions as an rms sum. 

"State-independence" means that the contributions are proportional to the 
wave function at the origin squared ^(O)) 2 or its first derivative 1^(0)^(0)1, 
i.e., to 5io/n 3 , with no other dependence on n. In particular, state-independent 
effects do not contribute to A(n) and D 2 \ considered in this paper. The non- 
leading terms, in contrast, depend on a state and involve sometimes substantial 
cancellations. They also may involve various numerical cancellations and their 
calculation can hardly improve the reliability of the estimation. However, their 
evaluation is necessary as a step in direction towards the completion of the 
calculation of the whole correction and to check whether the non-leading terms 
are of a reasonable value which sometimes may be enhanced 12 . We underline 
that it could be no rigorous proof for such an empirical estimation — it is 
barely a result of experience based on a high number of known logarithmic 
contributions. Such an approach is not much different from what experimen- 
talists do to estimate possible systematic errors. The estimation is always a 
guess, hopefully a plausible one, and definitely a deeper understanding of the 
problem is helpful. Our motivation is not only application of the factor of 1/5 
as a parameter, but also an observation that most functions in physics are 
quite smooth. The leading logarithmic term for some value of Za, which may 
be even much below 1/137 related to a non-physical value of Z <C 1, deter- 
mines the sign and the magnitude of the correction. In such a case even if the 
parameter approaches unity, the order of magnitude is often still determined 
by the leading asymptotics as long as the function is smooth enough. Most 
of such smooth functions seldom change their sign and that is also important 
for the motivation of our method. A weak point of this explanation is that we 
already deal with certain higher order corrections, i.e., with a function with its 
first few terms of expansion subtracted. Even if the function itself is a 'physi- 
cal' one with the expected smooth behavior, it is not quite clear how smooth 
it should be after the subtractions. Also, we know that some smooth functions 
change sign etc. and in this case the asymptotics may be not helpful. 

It seems that we have reached a certain numerical limit related to the log- 
arithmic contribution and the calculation of the non-logarithmic terms will 
be much more complicated than anything else done before. However, it is 
strongly needed to improve the reliability and the accuracy of the theoretical 



In particular, large coefficients were discovered for linear logarithmic term for the 
Lamb shift [67,68]. The coefficients are so large that the term linear in ln(l/(Za)) 
dominates over a cubic logarithmic term. The result has not yet been confirmed by 
an independent group. Neither origin of the large numerical value of the coefficient 
is understood, nor possible consequences for estimation of other higher-order terms. 
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predictions. 



13.3 Crucial 'soft QED' corrections 

Not only the crucial orders of magnitude offer a cross comparison of efficiency 
of different experiments. The modern bound state QED theory clearly recog- 
nizes two kinds of contributions: soft-photon contributions and hard-photon 
contributions. The latter are very similar to free QED, while the former essen- 
tially involve binding effects. There are two crucial soft-QED contributions. 

• One of crucial binding corrections is due to the higher order of the one-loop 
and two-loop self-energy. The one-loop self energy was recently success- 
fully calculated for most of important applications. The two-loop self-energy 
terms (see, e.g., Fig. 35) are partly known and uncertainty of this knowledge 
determines the accuracy of the Lamb shift calculations for hydrogen, helium 
ion, the medium- and high-Z Lamb shift and fine structure, while similar 
diagrams with a transverse photon contribute to an uncertainty of D 2 \ in 
the helium-3 ion. 



Fig. 35. One of diagrams for higher-order two-loop corrections (order of 
a 2 (Za) 6 mc 2 ) to the Lamb shift. The logarithmic part contains up to the cube 
of the logarithm. 

• The other important corrections are caused by recoil effects (see Figs. 36 
and 37) and they determine the uncertainty of the muonium HFS, the 
positronium energy levels (the hyperfine interval, the Is — 2s transition, 
and the fine structure) and a part of the uncertainty of the D 2 \ differ- 
ence for the helium-3 ion. The soft-photon nature of these contributions is 
clearly seen from the appearance of the bound logarithm \a(l/(Za)). Ac- 
tually this logarithm often (but not always) arises as ln(mc 2 / (E)) , where 
(E) is a characteristic atomic energy. The other possibility for this loga- 
rithm to appear is a logarithmic integration over coordinate or momentum 
space which leads to \n(mc{r) /h), where (r) is the characteristic size of the 
atomic state. Thus, the space-logarithm contributions contain a soft part 
with integrations over atomic momenta k ~ (Za)mc and hard part related 
to integration over k ~ mc. In energy- logarithm cases, two low-momentum 
regions (k ~ (Za) 2 mc and k ~ (Za)mc) can be distinguished. 

If one calculates separately the soft- or hard-momentum contributions, the 
logarithm of an effective energy or momentum appears as a logarithm of a 
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certain cut-off. That greatly simplifies evaluation of logarithmic corrections. 
Meantime, the appearance of such logarithm means that both the soft and 
hard contributions are divergent and only their sum have sense. 



Fig. 36. An example of the diagrams for higher-order radiative recoil corrections to 
the HFS interval (relative order of a(Za) 2 m/M). The logarithmic part contains up 
to the square of the logarithm. 



Fig. 37. An example of the diagrams for higher-order pure recoil corrections to the 
HFS interval (relative order of (Za) 3 m/M). The logarithmic part contains double 
logarithmic contributions. 



13.4 Crucial 'hard QED' corrections 

The hard-photon effects are very similar to effects of free QED, however, 
the areas of integrations over momentum space are different. Both the free 
QED and 'hard' contributions to the bound-state QED suffer from infrared 
divergencies. This problem is also solved in different ways. 

Most accurate free QED calculations which may be compared with experiment 
are related to the anomalous magnetic moments of an electron and a muon. 
The integration for them is performed over kind of an isotropic area in the 
Euclidean 4D space 13 , and the crucial level is a four-loop approximation. For 
the bound problems there are two other specific areas of integration. 

• For some problems it is sufficient to apply an external-field approximation 
to hard-photon corrections, in which the exchange photons carry zero en- 
ergy transfer (see, e.g., a diagram in Fig. 35). Even for recoil effects, the 

13 The physical expressions to calculate originate indeed from a description related 
to the Minkowsky space with the metrics (1, —1, —1, —1), where, e.g., k^k^ = A;q— k 2 . 
However, while calculating the diagrams it is often convenient to make the Wick 
rotation of the integration contour and to arrive at the Euclidian metric (1,1,1,1) 
(or rather (— 1, — 1, — 1, — 1)). 
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integration is essentially not covariant in the Euclidian 4D space. The high- 
est crucial orders are related to four-loop corrections for the external-field 
and three- loop corrections for recoil effects (see, e.g., diagrams in Figs. 36 
and 37). Still, for certain four- dimensional Euclidian variable one can ob- 
serve for characteristic momenta \ko\ <C \ki\. Note that in contrast to the 
anomalous magnetic moment, this is a calculation for two different particles 
and thus it has its own simplifications and difficulties. 
• The other specific situation for integration is related to the positronium 
annihilation when some photons are real and thus their momentum satisfy 
the condition k 2 = k% — k 2 = 0. Studies of the orthopositronium decay 
mode (see, e.g., Fig. 38) and a branching fraction for five-photon decay 
[253,194] allow us to check calculations of four-loop corrections in such 
a non-isotropic area of integration. The accuracy of determination of the 
branching of four-photon decay of parapositronium [190,191,192,193,194] is 
approaching a level where four-loop diagrams are important. 



Fig. 38. An example of the diagrams for two-loop radiative corrections to the 
three-photon annihilation of orthopositronium). 

The bound state problems supply us with an opportunity to check modern 
algorithms for four-loop calculations, competitive to the anomalous magnetic 
moment of an electron, and allow us to test them for different kinematic areas. 



14 What is next? 

Above, we have reviewed the present status on precision theory and experiment 
on light two-body atoms. In this section, we will briefly overview problems and 
possible experiments in the field. 

The Lamb shift was successfully measured in hydrogen, deuterium and an 
ion of He + . Experiments in neutral atoms have reached a certain level of 
accuracy, which will be difficult to improve very much without completely 
new approaches. However, a certain progress is possible and is expected from 
optical frequency measurements. 

Nevertheless, microwave experiments may be important. There remains a 
question about reliability of a microwave experiment by Sokolov and Yakovlev. 
Their results are reproducible within their experiment, which means that any 
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problem related to possible systematic effects could be, in principle, resolved. 
Their statistical uncertainty is competitive with the present optical experi- 
ments. Unfortunately, the experiment was stopped. We hope it may be con- 
tinued or reproduced. 

We also note that a microwave study of hydrogen circular states at MIT is only 
a project that can deliver a really independent value of the Rydberg constant, 
but, unfortunately, the evaluation of the data remains unfinished. 

In the case of the helium ion only microwave measurements have taken place. 
Ongoing optical experiments can offer essential progress. 

Tritium is a light neutral hydrogen-like atom whose the Lamb shift has not yet 
been measured accurately. We are interested in tritium not because of QED, 
but mostly because we hope to learn more about its nuclear structure and in 
particular on its charge radius. 

Uncertainty in determining the proton charge radius is the main obstacle in 
any improvement of the precision theory for the Lamb shift in hydrogen. The 
former scattering measurements were evaluated improperly, while their origi- 
nal data are not available. A new measurement, even with the same accuracy 
but with a transparent and reliable evaluation, would be very helpful. We can 
also learn the value of the proton charge radius by measuring the Lamb shift 
in muonic hydrogen; and we anticipate news from a related PSI experiment. 

The hyperfine interval of the ground state in hydrogen has been known for a 
while with a high experimental accuracy. Any theoretical progress is possible 
there only with a better understanding of the proton structure. In particular 
it is crucial to have a reliable description of the proton polarizability effects. 
This may be done with improvement of the experimental data on scatter- 
ing (muon-proton's or electron- proton's). We observed a certain interest by 
high-energy community to re-evaluate of various experimental data for the 
scattering (mainly due to determination of the proton charge radius). We 
hope for a reliable re-evaluation of existing data on inelastic scattering to be 
obtained soon. More data on spin-dependent inelastic muon-proton scattering 
are strongly needed. 

An advanced phenomenological description of the nucleon structure and nucleon- 
nucleon interaction, based on a constraint from experiment, is crucial not only 
for hyperfine splitting in hydrogen but also in few-nucleon light atoms, such 
as deuterium, tritium, helium-3 ion. Accuracy can be improved and we have 
observed some progress in this field. We hope that nuclear theory will help 
with prediction of the nuclear-structure contributions to the HFS interval for 
A = 3 nucleus. 

Successful comparison of the Is and 2s hyperfine intervals is not completely 
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sufficient. Accuracy of the experimental determination of the 2s HFS interval 
in hydrogen and deuterium is still not high enough, however, its threefold 
improvement after applying optical methods allows to hope for future progress. 
The 2s measurement in the 3 He + ion is still competitive with theory. However, 
it is not satisfactory when so important experimental value was never re- 
measured again. It seems that it is possible to go to somewhat higher Z and 
study Li ++ ions. 

The comparison of the Is and 2s HFS intervals is not the only way to get 
rid of the nuclear effects. Another important option is spectroscopy of pure 
leptonic atoms, such as muonium and positronium. Unfortunately, this kind of 
experiments was completely stopped some time ago. There are no running or 
planned muonium projects, while all present positronium projects aim only at 
the study of decay. We think that a measurement of the Is HFS interval, Is— 2s 
transitions in muonium and positronium and 2S — 2P fine structure intervals 
in positronium can deliver results quite competitive with theory. In the case 
of the hyperfine splitting in positronium the experiment is necessary also to 
resolve a discrepancy of approximately 3a which shows that our understanding 
of the situation is not completely satisfactory. 

Successful experiments designed to determine the g factor of a bound elec- 
tron in a hydrogen-like ion with a spinless nucleus eventually led to the most 
accurate determination of the electron-to-proton mass ratio (or, which is es- 
sentially equivalent, of the electron mass in unified atomic mass units). The 
experiment was done with carbon-12 and oxygen-16. Improvement in accu- 
racy is possible. It is also of crucial importance to present the data properly. 
A value g( 12 C 5+ )/g( 12 7+ ) should be presented by experimentalists by taking 
into account correlations between two measurements. The range of ions should 
be broadened to include ions from 10 Be 3+ to 40 Ca 19+ . Medium- Z lithium-like 
ions may be helpful at least to study systematic effects since they broaden a 
possible charge-to-mass ratio of ions and can be reasonably well understood 
theoretically. 

These experiments are more accurate than the MIT study of Breit-Rabi levels 
in hydrogen and deuterium which delivered the best values of /i p / /i e and /^/ A* e - 
However, these experiments were never reproduced and that is not satisfactory 
because of the great importance of an accurate value of fi p / /i e for QED tests 
and determination of fundamental constants. In principle a value of \i v j \i e can 
be and should be determined independently. 

One of possibilities can be a determination of the anomalous magnetic mo- 
ment 14 of the deuteron (similarly to a measurement of the anomalous mag- 

14 We follow the common notation which is somewhat misleading. What is called 
'anomalous magnetic moment' is rather 'anomalous g factor'. In this paragraph we 
are discussing the anomalous contribution to the magnetic moment. 
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netic moment of muon at the Brookhaven National Laboratory). That mea- 
surement was done for a value of a il -(ehjm v ) / fi p . A similar result for a deuteron 
would be very helpful since the deuteron anomalous magnetic moment is small 
(a d — —0.143) and its 'normal' value, related to the g factor gf^ = 2, i.e., 
gf^ ■ (eh/2m d ), in units of the nuclear magneton is well known (that is simply 
the mass ratio for a proton and a deuteron multiplied by two). Smallness of 
the magnetic moment provides an enhancement of accuracy in determination 
of the magnetic moments 

fXd_ = 2mp 2a d ■ (eh/2m d ) _ 
fx N m d lip Hn 

Comparing the deuteron and proton magnetic moments via NMR spectroscopy 
of HD molecules, one can find g p . 

Determination of the g factor of a proton (g p = h p /^n) is essentially the 
same as a determination of fi p /fi e , since the conversion factor [ji,n / He is known 
well enough. W. Quint suggested [254] to measure the proton g factor directly 
using technics developed for the experiment on hydrogen-like ions with spinless 
nuclei. 

The bound state QED theory is mainly satisfactory for hydrogen-like atoms 
for comparison with various present-day experiments. Taking into account 
future experiments, it is necessary to improve its accuracy calculating higher- 
order corrections. As one of the crucial problems, let us mention numerical 
and analytic calculations of one-loop contributions for the hyperfine difference 
D 2 i, two-loop corrections for the Lamb shift, D 2 \ and bound electron g factor 
and various fourth-order contributions to the positronium energy levels, the 
muonium Is hyperfine interval and D 2 \- 

This review is focussed on light atoms which offer various precision tests of 
bound state QED. "High precision" is understood as a high relative precision 
of various measurements and calculations and sometimes as a high sensitivity 
to certain high-order QED corrections. Certain study of medium- and high-Z 
ions can be also of interest for QED tests. However, it is crucially important 
to clarify uncertainty of the corrections due to the nuclear charge distribution. 

We note also that the high precision is not the only quantity of interest. One 
of the motivations to 'test' QED is a search for its violations; another is a test 
of various approaches for precision calculations and measurements. It is not 
necessary that the most precise tests are the most sensitive to any open ques- 
tions. The sensitivity is often quite selective and depends on a problem under 
question. Study of tests related to any particular problem needs a separate 
study and other experiments may be important. 
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15 Summary 



When they said, "Does it fit?" 
He replied, "Not a bit!" 

Edward Lear 

The situation concerning precision studies of simple atoms is fortunately oppo- 
site to that in the quoted verse by Lear. The theory and experiment agree well 
in general. There may be some discrepancy in a certain particular case, even a 
long standing one, but we have no reason to doubt quantum electrodynamics 
in its basic moments. However, the theory of simple atoms is a bound state 
QED theory and the crucial question is not whether QED proper is correct. 
The question is how successfully QED can be applied and has been applied 
to the bound states. The latter is not a simple issue and progress in QED 
calculations can be very fruitful for better understanding of the bound states 
in nuclear and particle physics. In early time of quantum physics theoretical 
progress established new theories and was of great significance. Now the most 
important part of physics of simple atoms lies in experiment. 

Fascinating experimental progress with high-precision measurements, by bridg- 
ing microwave, infrared, optical and in part ultraviolet domains of spectra, 
with dramatic cooling of atoms, with access to trapped single atoms and par- 
ticles, with producing muonic and exotic atoms and various highly charged 
ions offers us a broad range of applications. However, studying simple objects 
with predictable properties is still the most natural choice. That is how sim- 
ple atoms enter modern physics, most of progress of which is now rather a 
technological one and that is how simple atoms can serve for quite practical 
applications testing new experimental methods. Another application of sim- 
ple atoms of practical importance is the determination of certain fundamental 
constants. These constants have been more and more involved in the modern 
system of units and standards. 

Can we learn from simple atoms any really new physics? The answer for most 
experiments is rather negative. They have pursued other purposes. Simplicity 
of the atoms and reliability of the QED predictions make them a powerful tool 
to study effects beyond atomic physics, such as particle and nuclear properties. 
Precision study of two-body atoms delivers us information on structure of 
light nuclei (deuteron, triton, helion, and alpha particle), the structure of a 
proton, various properties of muons, pions and other particles. In general, we 
check consistency of QED theory and experiment and as long as we see the 
consistency, the constraints on new physics are related to the uncertainty of 
the experiment and theory. There are still some projects with simple atoms 
which are designed to look beyond the established physics: a search for exotic 
decay of positronium, a search for variations of the fundamental constants, a 
search for CPT violation studying properties of antihydrogen, a conversion of 
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muonium into antimuonium etc. 



Twenty years ago, when I joined the QED team at Mendeleev Institute for 
Metrology and started to work on theory of simple atoms, the accuracy of 
the experimental data for most QED tests was substantially higher than that 
of the theoretical predictions. Work of several theoretical groups from around 
the world during this twenty-year period has put theory of two-body atomic 
systems to a dominant position and we are now looking forward to obtaining 
new experimental results. Except the case of the 2s hyperfine interval in the 
3 He + ion, theory is more accurate than experiment. However, theory is still 
not perfect. For some quantities the gap between the theoretical and experi- 
mental accuracy is not big and we should not think that the theory is above 
any doubts. For example, starting my work on the review, I believed that the 
theory for the -D21 difference in helium-3 ion is more accurate than the exper- 
imental data, but a conservative re-estimation of the theoretical uncertainty 
[115,116] has reversed the situation. 

Several experiments are in progress or planned and substantial experimental 
progress is possible. Since the theoretical improvement needs some long-term 
programs, we have to start working on a further development of theory now. 
It will be indeed helpful for difficult theoretical calculations to be motivated 
by new experimental projects and new experimental results and I hope that 
both will follow. 

The obvious success in one-electron atoms turned attention of some theorists 
to three-body (such as helium, muonic helium, antiprotonic helium) and four- 
body systems and we hope that in some future it will be possible to state 
that theory of three- or maybe even four-body atoms has also superseded 
experiment. 
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n the principal quantum number 

I the orbital moment of an electron (in units of K) 

s the electron spin (in units of %) 

j the angular momentum of an electron (in units of ft) ( j = 1 + s) 

I the nuclear spin (in units of %) 

F the atomic angular momentum (F = j + I) 

m the mass of the orbiting particle (mainly, an electron) 

M the nuclear mass 

rriR the reduced mass of the orbiting particle (mj; = Mm/ (M + m)) 

Z the nuclear charge (in units of the proton charge) 

A the atomic mass number (i.e., M ss Am p ) 

i?oo the Rydberg constant (-Roo = a 2 m e c/2h) 

a the fine structure constant (a = e 2 /Aireohc) 

c the speed of light 

h the Planck constant 

h the reduced Planck constant (h = h/2it) 

e the proton (positron) charge 

fiB the Bohr magneton (fJ-B = eh/2m e ) 

fiN the nuclear magneton (fix = eh/2m p ) 

m e the mass of an electron* 

Hh the magnetic moment of a helion* (the nucleus of 3 He) 

a e the anomalous magnetic moment of an electron* (a e = (g e — 2)/2) 

Qp, the g factor of a muont (g^ = 2^/(e^/2m M )) 

Rn the nuclear charge radius* 

Re also the nuclear (rms) charge radius, when it is necessary to distinguish it from Rm 

Rm the nuclear (rms) magnetic-distribution radius 

Table A.l 

The most frequently used notation. 

* - similar for other particles and nuclei; * - similar for a muon; * - similar for an 
electron and a proton; * - similar for a proton, a deuteron etc. 
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A Notations 



Unless otherwise stated, we use SI units. However, all numerical results for 
energy E are presented in frequency units E/h. 

Energy levels in a single-electron atom are labelled as nlj or if the hyperfine 
effects are involved as nlj(F). The orbital momentum is presented by letters s 
for I — 0, p for for I — 1, d for I — 2, / for I — 3 etc. For few-electron systems 
(and positronium) the capital letters for angular momenta are used such as 
n K L j where k — 2S + 1. 



B Extract from the list of the recommended fundamental con- 
stants (CODATA, 2002) 

Here, we summarize the most recent CODATA recommended values of fun- 
damental constants related to QED. We need to emphasize that Table B.l 
contains results obtained in the 2002 15 adjustment of the fundamental con- 
stants [23], while through out the paper we compare original results with the 
older set of recommended values [21]. A reason for that is that a substantial 
part of results reviewed here was obtained after 1998 and was not included 
in the earlier adjustment [21] and thus comparing the results we clearly see 
the recent progress in the field. The recent adjustment [23] already included 
nearly all original results reviewed in this work. 



5 We remind that 1998 and 2002 are years of the deadline for collecting the data 
for the evaluation. 
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Table B.l 

New recommended values of the fundamental constants (CODATA, 2002, [23]) re- 
lated to atomic physics and QED. 

* The result for the Rydberg constant for the hydrogen atom Rh was not presented 
in [23] directly and it is derived from values for R^ and m e /m p . * We ignore in our 
paper (except of consideration of L>2i) the negative sign of the ratio of the magnetic 
moments. 
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